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ABSTRACT 

Among the 25 planetary systems detected np to now by gravitational mi- 
crolensing, there are two cases of a star with two planets, and two cases of a 
binary star with a planet. Other, yet nndetected types of triple lenses inclnde 
triple stars or stars with a planet with a moon. The analysis and interpretation 
of snch events is hindered by the lack of nnderstanding of essential characteristics 
of triple lenses, snch as their critical curves and caustics. We present here analyt¬ 
ical and numerical methods for mapping the critical-curve topology and caustic 
cusp number in the parameter space of n-point-mass lenses. We apply the meth¬ 
ods to the analysis of four symmetric triple-lens models, and obtain altogether 9 
different critical-curve topologies and 32 caustic structures. While these results 
include various generic types, they represent just a subset of all possible triple¬ 
lens critical curves and caustics. Using the analyzed models, we demonstrate 
interesting features of triple lenses that do not occur in two-point-mass lenses. 
We show an example of a lens that cannot be described by the Chang-Refsdal 
model in the wide limit. In the close limit we demonstrate unusual structures 
of primary and secondary caustic loops, and explain the conditions for their oc¬ 
currence. In the planetary limit we hnd that the presence of a planet may lead 
to a whole sequence of additional caustic metamorphoses. We show that a pair 
of planets may change the structure of the primary caustic even when placed far 
from their resonant position at the Einstein radius. 

Subject headings: gravitational lensing: micro — methods: analytical — plane¬ 
tary systems 
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1. Introduction 

In the past two decades gravitational microlensing surveys have been very successful, 
in particular as a tool for studying the stellar population toward the Galactic bulge. In a 
microlensing event, a star passing close to the line of sight to a background “source” star is 
detected by its gravitational lens effect, which temporarily amplihes the flux from the source 


e.g., Paczyhski 1996). The main advantage of the method is its sensitivity to low-mass 


objects, ranging from stellar down to planetary masses, with most of them too faint to be 
routinely detected by other means. 

In addition to single-star microlensing events, right from their start of operations the 


surveys have detected events with binary-star lenses (Udalski et ah 1994, Alard et ah 1995). 


The frequency of binary events is lower than the frequency of binary stars, since binaries with 
too close or too far components often mimic single-lens events. The microlensing sensitivity 
to low mass ratios dually led in 2003 to the first detection of microlensing by a star with a 
planet ( Bond et al.||2004[ ). By the time of writing, altogether 25 planetary systems had been 
detected by microlensing A majority involved a star with a single planet; nevertheless, 
four of them involved three-body systems. In two cases the lens was a star with two planets 


Gaudi et ah 

2008 

Han et al. 

2013) 

, while the other two involved a binary with a planet 

Gould et al.||2014 

Poleski et al.||201^ 

:). Other possible triple-lens systems that had not been 


detected yet include triple stars, or even lenses formed by a star with a planet with a moon. 

Turning to the underlying physics, in the case of a single lens the light from the source 
star is split into two images, which remain generally unresolved due to their small angular 
separation < 1 mas. The accompanying temporary amplihcation of the flux from the source 


typically produces a simple symmetrically peaked light curve (e.g., Paczyhski 1996). Lenses 
with multiple components produce a higher number of images and lead to a greater diversity 
of light curves, which peak anytime the source crosses or approaches the caustic of the lens. 
Any simple caustic crossing leads to the appearance or disappearance of a pair of unresolved 
images at positions dehning the critical curve of the lens. The caustic and the critical curve 
thus are key characteristics of the lens. They depend sensitively on the lens parameters: the 
masses and positions of the components. Understanding the range of possible critical-curve 
and caustic geometries is a prerequisite for successful analysis and interpretation of observed 
microlensing light curves. 

Microlensing by a two-component lens, such as a binary star or a star with a planet, 
is well described by the two-point-mass lens model. Such a lens has only two relevant 


Tt tp: / / exoplanet archive. ipac. caltech.edu 








































3 


parameters: the mass ratio, and the projected component separation. The model has been 


analyzed in detail by Schneider & Weiss (1986), Erdl & Schneider (1993), and Witt & 


Fetters (1993); its limiting cases were stndied by Dominik (1999). For any mass ratio the 


critical curve was shown to have three different topologies. In order of decreasing component 
separation these are: “wide” with two separate loops, “intermediate” (or resonant) with a 
single loop, and “close” with an outer plus two inner loops. Each regime has a corresponding 
caustic geometry with the same number of separate non-overlapping loops: wide has two 
four-cusped loops, intermediate a single six-cusped loop, and close has one four-cusped plus 
two three-cusped loops. A source positioned outside the caustic has three images, while a 
source positioned inside the caustic has five images. 

Microlensing by triple lenses can be described analogously by the three-point-mass lens 
model. There are five lens parameters: two relative masses, and three relative positions 
defining the two-dimensional configuration of the components in the plane of the sky. In 
comparison with the two-point-mass lens, the model has a number of qualitative differences. 
For example, varying the lens parameters may lead to a change in the cusp number of the 
caustic without any accompanying change in the topology of the critical curve. Such changes 


in the caustic structure occur via swallow-tail or butterfly metamorphoses (Schneider et ah 


1992). In addition, loops of the caustic may overlap, and individual loops may self-intersect. 


As a result, the caustic may have inner multiply nested regions. All caustics separate the 
outer four-image region from an inner six-image region. Only caustics with overlapping 
loops or self-intersections have additional eight-image regions, in case of double nesting even 
ten-image regions. 

However, unlike in the two-point-mass lens case, the full range of critical-curve topolo¬ 
gies and caustic structures of the three-point-mass lens has not been explored yet. Such a 
study would require a systematic mapping of the five-dimensional parameter space, detecting 
changes in the critical curve and caustic of the corresponding lens. Nevertheless, a range of 
published works have explored different specific regimes of triple lenses. The first studies 
explored binary lenses with an additional external shear (Grieger et al.||1989 Witt & Fetters 


1993), with caustics already displaying swallow tails and butterflies. Most numerous are the 


studies involving a stellar lens with two (or more) planets. Some of them demonstrate effects 
on the light curve, others for example place constraints on the presence of a second planet in 
observed single-planet events. Without claiming completeness, we refer here to the works of 


Gaudi et al. 

(1998) 

Bozza (1999; 

, Han et al. (2001), 

dan & Fark 

(2002) 

Rattenbury et al. 

(2002), Han 

(2005) 

Kubas et al. 

2008), Asada (2009) 

, Ryu et al. 

2011|, 

Song et al. (2014), 


and Zhu et al. (2014) 


Lensing by a binary star with a planet has been explored less frequently ([Bennett et 
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al. 

1999 

Lee et al. 

2008 

Han 

2008a 


2010|. In view of the two recently 


detected systems, this is bonnd to change. The close and far limits of triple-star lensing were 
investigated by Bozza (2000a|[b). Finally, lensing by a star with a planet with a moon was 


stndied by Han & Han (2002), Gandi et ah (2003), Han (2008b), and Liebig & Wambsganss 


(2010), with prospects for detecting snch systems remaining open. 


Here we set ont to systematically stndy the critical cnrves and canstics of triple lenses. 


Based on the work of Erdl & Schneider (1993) and Witt & Betters (1993), we develop methods 


for efficient mapping of critical-cnrve topologies and canstic geometries in the parameter 
space of n-point-mass lenses. We then apply the methods to the analysis of simple triple¬ 


lens models. Initial steps of this research appeared in Danek (2010), Danek & Heyrovsky 


(2011), and Danek & Heyrovsky (2014). 


We start in §|^by introdncing the basic concepts of n-point-mass lensing. In particniar, 
we concentrate on the Jacobian and its properties, snch as the eqnivalence of its contonrs 


with critical curves of re-scaled lens configurations (Danek & Heyrovsky 2015). Analytical 


and numerical methods for mapping critical-curve topologies in the lens parameter space 
are introduced in §i In§i we discuss the caustic and its metamorphoses, and show how 


to track changes in cusp number using the cusp and morph curves (Danek & Heyrovsky 


2015). We apply the methods to triple lenses in § [^ starting with a brief overview of their 
properties in § |5.1[ In § 5^ - § |5.5| we include a full analysis of four symmetric two-parameter 
triple-lens models, with an overview of the found critical curves and caustics in § 5.6 We 
end by summarizing the main results and highlights in § [^ 


2. The n-point-mass lens and its Jacobian 


Galactic gravitational microlensing events can be described using a simple n-point-mass 
lens model, consisting of n components (stars, planets) in a single lens plane with no external 
shear and no convergence due to continuous matter. Following Witt (1990) we describe 
positions in the plane of the sky as points in the complex plane, with separations measured 
in units of the Einstein radius corresponding to the total mass of the lens. 


The relation between the position of a background source ( and the position of its 
image formed by the lens is expressed by the lens equation 


n 



( 1 ) 


where Zj and /ij are the positions and fractional masses, respectively, of the individual 
components of the lens, and bars over variables denote their complex conjugation. The 
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complex plane of the image positions 2 ; is called the image plane, and the complex plane of 
the source positions ( is called the source plane. The fractional masses are normalized to 
the total mass of the lens, so that hi = 1- 

The positions of individual components zj generally change with time as they orbit 
around the center of mass of the lens. In this work we study the lensing properties of 
a “snapshot” n-body conhguration at a given instant. The case when the change of lens 
parameters is non-negligible on the timescale of the microlensing event can be described by 
a corresponding sequence of such snapshot conhgurations. 


The n-point-mass lens with n > 1 components produces between n + 1 and 5 (n — 1) 
images (in steps of two) of any point in the source plane (Rhie||2003 Khavinson & Neumann 


2006). In the image plane images appear and disappear in pairs along the critical curve of 
the lens. The critical curve can be expressed as the set of all points Zcc for which the sum 




hi 


ZjY 


= e 


-2ic 


( 2 ) 


lies on the unit circle (Witt 1990| . Here 0 is a phase parameter spanning the interval 
0 G [0,7r). In the source plane the image number changes when the source crosses the 
caustic Cc of the lens. The caustic is obtained by tracing critical-curve points back to the 
source plane using equation ([^, i.e.. 


Cc ^cc ^ ^ 


hi 


( 3 ) 


i=i 


In mathematical terms, the critical curve is the set of points in the image plane with 
zero Jacobian of the lens equation. The Jacobian 


det J (z) = 1 


n 

E 

i=i 


hi 


( 4 ) 


is discussed in detail in Danek & Heyrovsky (2015). Here we only summarize its properties 


important for the following analyses. As seen from equation (|^, the Jacobian is a real 
function dehned over the complex plane, running from —00 at the positions of all components 
z = Zj to 1 at complex inhnity and at the positions of all Jacobian maxima. These can be 
found as the roots of the polynomial obtained from 


E hj 

(z — Zj)^ 
i=i ^ 


0 . 


( 5 ) 
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The degree of the corresponding polynomial indicates the Jacobian has up to 2 n —2 different 
maxima. The number may be lower if there are any degenerate roots; these correspond to 
higher-order maxima. A doubly degenerate root corresponds to a double maximum, a root 
with degeneracy 3 corresponds to a triple maximum, etc. 


The saddle points of the Jacobian can be found similarly among the roots of the poly¬ 
nomial obtained from 

n 

V , .3 ^0. (6) 

The corresponding polynomial has up to 3 n —3 different roots. First we need to sort out any 
common roots of equations ([^ and (|^, those correspond to higher-order maxima instead of 
saddles. All remaining roots are Jacobian saddle points. The number of different saddles 
may be reduced further if there are any degenerate roots; these identify higher-order saddles. 
A non-degenerate root corresponds to a simple saddle, a doubly degenerate root corresponds 
to a monkey saddle, etc. 


Studying the contours of the Jacobian, Danek & Heyrovsky (2015) pointed out a re¬ 
markable correspondence. While the zero-Jacobian contour is the critical curve Zcc of the 
lens, any other Jacobian contour zx with det J (2;a) = A is a re-scaled critical curve of a lens 
with the same components in re-scaled positions. Denoting by Zcc{fij,Zj) and zx^fij^Zj) the 
critical curve and det J = X contour, respectively, of a lens with masses /ij and positions Zj, 
we can express this correspondence by 


^cc(h'j)\/l A) / \/l A. (7) 

For A = 0 we get the critical curve of the original conhguration. Close to the original 
positions, the A —)■ — oo contours are shrunk versions of the wide-limit critical curves. At the 
other limit, the highest-Jacobian A —>■ 1 contours are expanded versions of the close-limit 
critical curves. A single Jacobian contour plot thus yields the critical curves for all scalings 
of the given lens configuration, from the close to the wide limit. 


3. Critical-curve topology regions in parameter space and their boundaries 

3.1. Critical-curve topology and its changes 

Here we summarize general properties of n-point-mass lens critical curves following 

. Varying the phase parameter 0 in the critical curve equation ([^, 
we obtain 2 n continuous line segments, which may connect to form 1 < Nioops <2n closed 


Danek & Heyrovsky (2015 
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loop^ of the critical curve (Witt 1990). Individual loops may lie in separate regions of 
the image plane, or they may lie nested inside other loops. The total number and mutual 
position of loops dehne the topology of the critical curve. 


In the wide limit, all lenses have a critical curve with Nloops = n separate loops, corre¬ 
sponding to n Einstein rings of the individual components. In the close limit, the critical 
curve has Nioops = 1 + N^ax loops, where N^ax is the number of different Jacobian max¬ 
ima. One loop corresponds to the Einstein ring of the total mass, plus there is a small loop 
around each maximum of the Jacobian. In view of the discussion following equation ([^, 
Nmax < 2 n — 2, SO that in the close limit Nioops < 2 n — 1. The equality holds if the Jacobian 
has only non-degenerate (simple) maxima. Any degenerate (higher-order) maximum reduces 
the number of loops. For example, the two-point-mass lens always has two different maxima, 
thus its critical curve always has three loops in the close limit. 


When varying the parameters of the lens such as its scale, the topology may undergo 
changes when individual loops merge or split. This occurs at Jacobian saddle points when the 
critical curve passes through them (e.g., Erdl & Schneider 1993). Two loops come into contact 
at a simple saddle, three at a monkey saddle, more loops at gradually higher-order saddles. 
As shown by equation the Jacobian may have up to 3n — 3 different saddle points, 
with the highest number occurring when there are no higher-order saddles and no higher- 
order maxima. The number of different Jacobian contours passing through the set of saddles 
identihes the total number of changes in critical-curve topology encountered when varying 
the scale of the lens from the wide to the close limit. This can be seen as a consequence of 
the Jacobian-contour / critical-curve correspondence expressed by equation ([^. Therefore, 
the critical curve of an n-point-mass lens may undergo no more than 3n — 3 changes in 
topology between the wide and close limits. 


In the case of the two-point-mass lens, there are always three simple saddles. One lies 
on the axis between the components, while an off-axis pair of saddles lies on a different 
Jacobian contour. The two distinct saddle contours imply that the two-point-mass lens 
critical curve always undergoes two changes when varying the component separation s, and 
thus has exactly three topologies. The wide topology has two separate loops, the intermediate 
topology has a single merged loop, and the close topology has an outer loop plus two small 
inner loops around the Jacobian maxima. The topology sequence is independent of the 
second lens parameter, the mass ratio of the lens components. 


Proceeding to lenses with more than two components, we note that the shape of the 


^We note that at least for n < 3 the sharp upper bound is 2 n — 1, and even for triple lenses we have 
found no more than 5 loops so far. 








critical curve depends on 3n — 4 lens parameters (e.g., Danek & Heyrovsky 2015). Follow¬ 


ing the preceding discussion, boundaries between regions in parameter space with different 
critical-curve topology can be found by identifying parameter combinations, for which the 


critical curve passes through a saddle point of the Jacobian (e.g., Erdl & Schneider 1993). 
The search for these boundaries is thus mathematically reduced to finding the conditions for 
the occurrence of a common solution of equations ([^ and (|^ . 

The usual analytical approach described in the following § |3.2 is based on rewriting both 
equations in polynomial form and computing their resultant. This step is then followed by 
a second resultant constructed from the Erst. However, this approach often yields unwieldy 
expressions. In addition, one has to check the results for spurious solutions. These do not 
occur in the two-point-mass lens, but they do appear in all the triple-lens models studied 
further below. 


As an alternative, we present in § 3.3 and § 3.4 two efficient numerical methods for 


mapping the boundaries. These can be used in models with at least one scale-defining 
parameter, such as the component separation s in the two-point-mass lens. In the first 
method we hnd the roots of the hrst resultant condition, while in the second we utilize the 


Jacobian scaling properties described by Danek & Heyrovsky (2015). Both methods are free 
of spurious solutions. 


3.2. Analytical boundaries computed by resultant method 


The method described here was pioneered in the context of critical-curve topology map¬ 
ping by Erdl & Schneider (1993) and Witt & Betters (1993). Multiplying the saddle-point 
equation ([^ by the product of its denominators yields 


Psadd{^} ^ ^ Pj (•^ ^k) 0 , 

j=l k=l,k^j 


( 8 ) 


a polynomial of degree 3 n — 3. In a similar manner we convert the critical-curve equation (|^ 
to 

n n n 

Pcritiz) = JJ (2: - Zkf JJ {z- Zkf = 0 , (9) 

k=l j=l k=l,k^j 

a polynomial of degree 2n for any value of the parameter 0. The analytical condition for 
the existence of a common root of Psadd{z) and Pcrit{z) is 


^GSzi^Psaddy Pcrit) 0 ) 


( 10 ) 
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where the resultant IleSz{f,g) of two polynomials f,gissi function of their coefficients. It 
may be computed by evaluating the determinant of the Sylvester or Bezout matrices, as 
described in Appendix 


The expression obtained from equation (10) is a polynomial in terms of 
denote tn = we can write the result as 


If we 


Presiw) = = 0 , 

j=0 


( 11 ) 


where the degree of the polynomial m < 3n — 3. The boundary condition is now equivalent 
to the condition for p^es to have a root on the unit circle. 


In order to obtain the condition purely in terms of the lens parameters, [Witt fc Fetters 


(1993) suggested the following approach for the two-point-mass lens. For a root along the 


unit circle = w \ so that if we take the complex conjugate of equation (11) and multiply 
it by tn™', we get another polynomial equation 


Pconjig^^ ^ ^ Otm—j ^ 0 

j=0 


( 12 ) 


Along the boundary in parameter space, polynomials pres and Pconj must have a common 
root, thus 

ReS^o (pres, Pconj) 0 (13) 


yields the sought condition in terms of lens parameters. We point out that equation (13) 


presents a single constraint in parameter space. Hence for two-parameter models (such as 
the two-point-mass lens or the triple-lens models described in § 5.2 - § 5.5) it generally 


describes a set of curves, for three-parameter models a set of surfaces, etc. 


We illustrate the approach here on the case of the two-point-mass lens (Erdl & Schneider 


1993 Witt & Fetters 1993), parameterized by the fractional mass of one component p G (0,1) 
and the separation between the components s > 0. If we align the lens with the real axis 
and place the components symmetrically about the origin, their positions and masses are 
{^ 1 , 2 : 2 } = {—s/2, s/2}, and {pi,p 2 } = (h, 1 — p}- The polynomial equations for Psadd and 

Pcrit 


3 3 1 

Psaddiz) = + -s{l- 2ij)z‘^ + -s^z + -s^(l - 2p) = 0 


and 


Pcrit(z) = z"^- -(s^ -F 


s(l-2ju) e^'’^z + 


1 

16 


= 0 , 


(14) 


(15) 


Using them in equation (10) leads to 

ReSziPsadd,Pcrit) = ” p)^ - 3s^(l - 9p + 9p^) -1- 3s‘^e^‘'^ - s®] =0. (16) 
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The factors in front of the square brackets are non-zero for any genuine two-point-mass lens. 
If we set = w, the term in the brackets yields the polynomial equation 


Pres{w) = — 3s^(l — 9fi + — s® = 0 . 


Following equation (12) we construct 


Pconj{w) = —s^w^ + — 3s^(l — 9^ + 9fi^)w -|- 1 = 0 


(17) 


(18) 


The resultant obtained from equation (13) can be factorized as follows: 


[1 - 3 s2(l - 9/i + V) + 3 + 3 _ 9;^ + 9^2^ + 3 + s®] 

X [1 - 3 / + 3 /(I - 9/i + 9/1^) - ^ = 0 . (19) 

At least one of these three square brackets thus has to be equal to zero. 


The first bracket in equation (19) is equal to Presi^), with w = 1 corresponding to 0 = 0. 
Therefore, it must include any transitions occurring on the critical curve along the real axis, 
in this case the passage of the critical curve through the central saddle point. Taken as a 
polynomial in s, the hrst bracket has a single real positive root (Erdl & Schneider 1993) 

-1 3/2 


— 


~ p 


( 20 ) 


which is the boundary between the wide and intermediate topologies. 


The second bracket is equal to —pres(—1), with w = —1 corresponding to 0 = 7r/2. 
There is only a single zero point of this expression in parameter space, [p, s] = [0.5, \/(h5], 
which corresponds to the intermediate - close splitting along the imaginary axis of the critical 
curve of an equal-mass lens. 

The third bracket corresponds to topology transitions that occur at any other values of 
0 on the critical curve. Here they describe the passage of the critical curve through the pair 
of saddle points off the real axis. Taken as a polynomial in s, the third bracket has a single 
real positive root (Erdl & Schneider 1993; Rhie & Bennett 1999) 


Sn — 




- 3/4 


= 

' 5/1 


( 21 ) 


which is the boundary between the intermediate and close topologies. This boundary in fact 
passes through the single zero point of the second bracket in equation (19) at /i = 0.5 as 
well. The two curves given by equations (20) and (21) thus fully describe the division of the 
two-point-mass lens parameter space according to critical-curve topology. 
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Applying the same procedure to two-parameter models of triple lenses leads to the 
following problem. The method yields not only all the sought boundaries, but also curves 
corresponding to no change in critical-curve topology. These spurious results are additional 
solutions of equation ( [I^ , which correspond to common roots of Pres{w) and Pconj{w) lying 
off the unit circle. In Appendix we describe how this occurs and why it does not occur in 
the two-point-mass lens. All curves corresponding to specific values of w that can be factored 


out of the final equation (13), such as Pres(l) and Pres(~l) in the case of equation (19), are 
genuine boundaries. Parameters satisfying the remaining parts of equation (13) should be 


verihed either by computing sample critical curves, or by confirming the existence of a root 
of pres{w) on the unit circle. 

Alternatively, instead of computing the second resultant, one may solve the polynomial 
from equation ( [TT| directly by re-substituting w = cos 2(p + i sin 2 0. By taking the real and 
imaginary parts of pres separately, equation 0 may be treated as a set of two real equations. 
In case all the coefficients aj are real, such as in the two-point-mass lens or all the triple-lens 

the imaginary part of equation ([IT| instantly yields the 


models described in 5 5.2 


5.5 


0 = 0 and 0 = 7r/2 solutions. The corresponding real parts then yield the Pres(l) = 0 and 
Pres{—^) = 0 boundaries, respectively. The other boundaries can be found by substituting 
cos 2 0 = ^ and solving the two obtained real polynomial equations in Finding all roots 
.^0 £ [“Ij 1] of one of the polynomials and substituting them into the other polynomial yields 
the remaining boundary conditions. 


3.3. Numerical boundaries computed using the first resultant 

For any particular n—point-mass lens model the lens masses pj and positions Zj depend 
on a given set of parameters. Mapping the boundaries in parameter space numerically can be 
a tedious task. If we can compute analytically at least the first resultant and the polynomial 
Pres given by equation 10. the task is reduced to keeping track of the roots of Pres as 
one scans through parameter space. Positions at which the absolute value of any of the 
roots crosses unity correspond to points on the boundary. However, even this can be a very 
time-consuming exercise. 

In case at least one of the parameters of the lens model defines an angular scale in the 
plane of the sky (such as the separation s in the two-point-mass lens), it is possible to reduce 
the dimension of the space to be scanned by one. At the same time this approach directly 
yields all boundary points for each reduced parameter combination. In case there are several 
angular-scale parameters, we select one of them (further denoted by s) and instead of the 
others we use their ratio to s. We can now assume that all masses are independent of s 
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and all lens positions are directly proportional to s. All the models discussed in this paper 
satisfy these requirements. 

If we divide equation ([^ by we obtain a polynomial equation of exactly the same 

form with z' = zjs and re-scaled lens positions z'y. = z^/s, 

n n n n 

p'saddiz') = Pj n = n = 0 • (22) 

j=l k=l,k^j j=l k=l,k^j 

The equation is thus explicitly independent of s. If we similarly divide equation ([^ by 
we obtain 

n n n 

Pcritiz) = Y[iz' - z'kf - s-2 Y Pj n = 0 , (23) 

k=l j=l k=l, k^j 

a polynomial of exactly the same form as Pcrit{z) but with s~^ instead of Setting 
the resultant of p'saddiz') and p'crit{z') equal to zero, we obtain a polynomial in terms of 

yjf — g2l(/l ^-2^ 

m 

Pres (w') = a'- = 0 , (24) 

j=0 

in which all coefficients a) are independent of s, and the degree of the polynomial m < 3 n—3. 
Finding the boundaries now requires a scan of the reduced parameter space (skipping s), 
with the boundaries directly given by all roots of p'^^^iw'). Any such root yields the 
value of s at the boundary, s = the same time it also yields the position along 

the critical curve, 0 = /2, at which the critical curve splits. 

Furthermore, if all coefficients a'- are real (such as in all the models discussed in this 
paper), all roots of p'resi'^') either are real or they occur in complex-conjugate pairs. A 
complex-conjugate pair of roots yields the same boundary value of s, corresponding to a 
simultaneous split at complex-conjugate positions along the critical curve. It is thus useful 
to keep track of the discriminant of p'^^g{w') in the reduced parameter space. Its change 
of sign corresponds to a transition from a pair of real roots (two boundaries) to a pair of 
complex-conjugate roots (one boundary). Zero discriminant corresponds to a multiple root 
of p(.g^(tc'), describing a single boundary point and a single 0 value along the critical curve. 

For illustration, for the two-point-mass lens we obtain 

p'^^^{w') = w'^ - 3(1 - 9/i + + 3 w' - 1 = 0 , (25) 


and its discriminant 


A{p'„,) = -3‘>i,^(i-rr(i-2ry<0. 


( 26 ) 
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The only non-degenerate configuration with zero discriminant has fi = 1/2, for which there 
is a double real root plus another real root: = {—2, —2, 1/4}. These correspond 

to boundary points s = {y/l/2, ■\/l/2, 2}. For all other /i G (0, 1) the discriminant is 
negative. Thus, has one real root and a complex-conjugate pair of roots, resulting 

in two boundary values of s for any fi. 


3.4. Numerical boundaries computed using Jacobian-contonr correspondence 


We return here to the method mentioned briefly in § 3.1 Instead of computing the re¬ 


sultant, we can map the boundaries in parameter space utilizing the correspondence between 


Jacobian contours and critical curves of re-scaled lens configurations (Danek & Heyrovsky 


2015), described here at the end of § 2 Just like in the previous case, even this method 


requires having a single angular-scale defining parameter s, with all lens positions directly 
proportional to it, and all masses independent of it. 


For a given set of the remaining parameters plus an arbitrary non-zero scale parameter 
So we compute the positions of the Jacobian saddle points Zsadd by finding all roots of Psaddiz) 
from equation (|^. For each root we evaluate det J (zgadd) from equation (|^, which identifies 
the contour passing through the saddle point. Now the Jacobian-contonr / critical-curve 
correspondence tells us, that in order for the critical curve to pass through the saddle point, 
we need to re-scale the lens positions by -^1 — det J {zsadd), as shown in equation Q. Hence, 
the boundary value of the parameter s is So-^l — det J (zgadd) ■ Recalling the polynomial 
degree of Psadd(z), this approach yields up to 3 n — 3 different boundaries for any given point 
of the reduced parameter space. 


For illustration we refer to Figure 1 from Danek & Heyrovsky (2015), which shows a 
Jacobian contour plot of a binary lens with p = 0.8 and Sq = 1- The real saddle point lies on 
the det J ~ —11.0 contour, and the complex-conjugate saddle-point pair lies on the det J ~ 
0.711 contour. The corresponding boundary separations are thus s = ^1 J- 11 ~ 1.861 and 
s = ^1 — 0.711 0.733, respectively. 


4. Mapping the number of cusps of the caustic in parameter space 

The caustic of the lens, defined here by equation ([^, consists of the same number of 
closed loops as the critical curve. We refer to Danek & Heyrovsky (2015) for an overall 


discussion of the properties of the n-point-mass lens caustic and its cusps. Here we are 
interested in tracking changes in the number of cusps in the parameter space of the lens. 
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Since individual loops or parts of the caustic may overlap in lenses with three or more 
components, we specify that we count the number of cusps encountered as we trace along 
each loop of the caustic. This means that a cusp superimposed over a fold is still a cusp, or 
that two or more cusps along different parts of the caustic that happen to lie at the same 
position in the source plane are still counted as two or more cusps. 


The cusp number changes in caustic metamorphoses. The simplest of them is the beak- 
to-beak metamorphosis, in which two tangent folds reconnect and form two facing cusps 


see Figure 5 in Danek & Heyrovsky 2015). In the image plane this always corresponds 


to a change in critical-curve topology at a saddle point. The corresponding boundaries in 
parameter space are thus identical to those studied in § |^ 

While the beak-to-beak metamorphosis already occurs in two-point-mass lenses, two 
additional types occur in triple lenses: the swallow-tail and butterfly metamorphoses. They 


are discussed in more detail and illustrated in Danek & Heyrovsky (2015). Both occur on the 


caustic without any change in topology or other significant effect on the critical curve. The 
lens-parameter conditions for the occurrence of these additional metamorphoses form new 
boundaries in parameter space, which subdivide the critical-curve topology regions discussed 


m 


For the following analysis we will use the cusp curve and morph curve introduced by 


Danek & Heyrovsky (2015). The parametric form of the cusp curve is 


_P 

(y. - 
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= A 
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(27) 


where A > 0 is a real parameter. The intersections of the cusp curve with the critical curve 
in the image plane identify the positions of cusp images along the critical curve. Since the 
curve is scale-invariant, its intersections with other Jacobian contours identify the positions 
of cusp images for all arbitrarily shrunk or expanded lens conhgurations. 


Another scale-invariant curve in the image plane is the morph curve, defined in para¬ 
metric form by 
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where F is a real parameter. Its intersections with the cusp curve away from the lens 
positions and higher-order maxima identify the positions of metamorphosis-point images. 
Intersections at saddles correspond to beak-to-beaks, simple intersections away from special 
points correspond to swallow tails, and intersections at two-branch self-intersection points 
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of the cusp-curve correspond to butterflies if the critical curve is tangent to either of the 


branches (Danek & Heyrovsky 2015). 


In order to hnd conditions for the occurrence of swallow-tail and butterfly metamor¬ 
phoses, we hrst locate all intersections of the cusp curve and morph curve, and leave out the 
lens positions, saddle points, and higher-order maxima. We use the parametric polynomial 


forms of the curves obtained by multiplying equations (27) and (28) by all their denomina¬ 


tors. We then employ a method analogous to the one used in § 3.2 for hnding the topology 


boundaries in parameter space. We compute the resultant of the two polynomials and ex¬ 
press the result as a polynomial in the morph-curve parameter F. Since we seek solutions 
with r real, we compute a second resultant of the obtained polynomial with its complex 
conjugate. 

The result can be factorized to obtain several conditions involving the non-scaling pa¬ 
rameters of the lens, each as a polynomial in the cusp-curve parameter A. For any combina¬ 


tion of the non-scaling parameters we look for real non-negative roots A. Equation (27) then 


yields a set of z’s on the cusp curve: those that satisfy equation (28) and he on the morph 


curve are the sought intersections. The hnal check is important, since the second resultant 
condition may introduce additional spurious solutions. In case the second resultant cannot 
be simply factorized or generally is too unwieldy, one can hnd the intersections numerically 
by gradually varying A and tracing the branches of the cusp curve, checking each point 2 ; 


with the morph-curve equation (28) 


We then compute the Jacobian at the metamorphosis points to identify the critical curve 
passing through each of them. The corresponding scaling factor \/Y 

The 


det J determines 

the separation of the lens components, just like in the saddle-point case in 


3.4 


obtained parameter combinations form additional curves in parameter space, leading to its 
hnal division into subregions with diherent combinations of cusps on the loops of the caustic. 


5. Triple-lens models 
5.1. General properties 

The general triple lens is described by hve parameters: typically two relative masses 
and three spatial-conhguration-dehning parameters. Among these three it is advantageous 
to select a single scale parameter, such as the separation of two of the components, or the 
perimeter of the triangle connecting the components. The remaining two parameters can be 
for example angles or relative lengths. With such a choice one may utilize the methods from 
§ and § 1^ for studying critical curves and caustics based on the scaling properties of the 
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Jacobian. 


The Jacobian surface has three poles at the lens positions, hence the critical curve has 
three loops in the wide limit. There are up to four different Jacobian maxima, therefore 
in the close limit the critical curve has up to hve loops: an Einstein ring corresponding to 
the total mass, plus up to four small loops around the maxima. The surface has up to six 
different saddle points. As a consequence, the critical curve of any triple-lens conhguration 
may undergo up to six topology changes when varying the scale parameter from the close to 
the wide limit. 


The critical-curve polynomial Pcrui^) from equation ([^ is of sixth degree, implying that 
the critical curve and caustic may each have no more than six loops. We point out here in 
advance that the models discussed further in 


5.2 


5.5 have only up to hve loops. 


A compact expression for Pcrit{z) can be obtained in terms of symmetric moments of 
the lens conhguration. In general, hve such moments are sufficient: 


k=l 


for j = 1... 3, and the mass-weighted 


dj — pk • 


k=l 


(29) 


(30) 


for j = 1,2. With these dehnitions di is the center of mass of the lens, and Ci/3 is its 
geometric center. In equal-mass models {pk = 1/3) the two centers coincide, so that di = 
ci/3, and d 2 = 2c2/3. 

If we set the origin of the complex plane at the geometric center (ci = 0), we arrive 
at a simpler form of equation (|^ than if we used the center of mass. The critical-curve 
polynomial equation can then be written using the four other complex moments: 


Pcrit{z) = - (2c2 + -2 (cs -h z'^ + [c\ + {2 C 2 - 3^2) 

+ 2c3{c2 - e^''^) z + cl- {cl + Cs di - C 2 ^ 2 ) = 0 .(31) 


In a similar manner we may re-write the Jacobian saddle-point polynomial equation (|^ 
in terms of the same four moments: 

Psaddiz) = J- 3 diZ^ — 3 (c2 — 2 ^2) Z^ J- (7 C3 J- C2 di) J- 3 (c2 J- 2 C3 (Ji — C2 ^2) z"^ 

— 3 Cs{c2 — ^ 2 ) ^ — cl{c2 — (^ 2 ) + cs{cs — C 2 di) = 0 . (32) 
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For any specific triple-lens model it is snfiicient to express the moments in terms of the model 
parameters and use them in equations (31) and (32). The transitions between different 


critical-curve topologies can then be mapped using the methods described in § 3.2 - § 3.4 


|5.5| we illustrate the nature of triple lenses by performing the 
analysis described in § |^ and § |^ on four simple two-parameter models. We chose the most 
symmetric models with a single scale parameter, sketched in Figure Two of the models 
are linear (symmetric with variable central mass, general asymmetric with equal masses), 
and two are triangular (equilateral with variable mass in one vertex, isosceles with variable 
vertex angle and equal masses). Each of the models, seen as 2D cuts through the full 5D 
triple-lens parameter space, intersects at least one other model. For example, the isosceles 
model includes equal-mass equilateral and equal-mass symmetric linear conhgurations. As 
a combination, these models form a reference set for further studies of the triple lens. 


In the following § 5.2 


For all the models we present parameter-space maps of critical-curve topologies and 
cusp numbers, together with galleries of all different topologies. In the accompanying text 
we pay particular attention to interpreting the wealth of information contained in the fig¬ 
ures. We point out any particularly interesting features of the critical curves and caustics. 
For better orientation, the text for each model is structured under the following headings: 
Model description, Jacobian surface character. Topology boundaries. Critical-curve topolo¬ 
gies, Caustic structure, Jacobian contour plots. Planetary limits. Close limit, and Wide limit. 
For the hrst model, the linear symmetric conhguration in § |5.2[ we additionally include the 


equations of the topology boundaries obtained from the resultant method, a gallery of all 
transition topologies with merging critical-curve loops, the polynomial expressions for the 
cusp and morph curves, and a gallery of caustics from all subregions of the parameter space 
subdivided by swallow-tail and butterfly boundaries. 


5.2. LS Model: Linear Symmetric Configuration 

Model description — in the LS model, sketched in the first panel of Figure [T| the three lens 
components he equidistantly along a line, spaced by separation s. The second parameter 
is the fractional mass pi of the central component. Placing the origin at the geometric 
center and the lenses along the real axis, the positions and masses of the components are 
{^ 1 , 2 : 2 , 2 ^ 3 } = {-s,0,s} and {/ii,/r 2 , p-s} = {(1 - Ai)/2,p, (1 - /^)/2}, with parameter ranges 
p G (0,1) and s G (0, cxd). 

The p = 0 limit describes an equal-mass binary with components 1 and 3 separated 
by 2 s. The p = 1 case corresponds to a single central lens represented by component 2. 
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The three-equal-masses case with /i = 1/3 is identical to special cases of two other models 
discussed further below: the p = 1/2 conhgurations of the LA model in 
conhgurations of the TI model in 


5.3, and the 9 = ti 


5.5 


Jacobian surface character — for fi ^ 1/9,1/5: four simple maxima -|- six simple saddles; 
for /i = 1/9: two double maxima -|- four simple saddles; for /i = 1/5: four simple maxima 
-|- two simple saddles -|- two monkey saddles. 


For fi < 1/9 all maxima lie in complex-conjugate pairs along the imaginary axis; for 
fi > 1/9 there are two complex-conjugate pairs of simple maxima symmetrically displaced 
from the imaginary axis. 


Two simple saddles lie along the real axis between neighboring lens components. The 
remaining saddles for fi < 1/5 he in complex-conjugate pairs along the imaginary axis; for 
fi > 1/5 they he in complex-conjugate pairs symmetrically displaced from the imaginary 
axis. 


Topology boundaries — due to the high degree of symmetry of the LS model, two of the 
moments defined by equations (29) and (30) are zero, c^ = di = 0. This leaves only two 
non-zero moments, C 2 = and d 2 = {1 — fi) s^. The critical-curve polynomial is reduced to 


Pcrit{z) = — {2s‘^ + z‘^ + [s"^ + {3fi — 1) s'^z'^ — fis‘^ = 0 , 


(33) 


and the saddle-point polynomial to 

Psaddii) ^ + 3{1 -2fl) + 3p s'^z^ - ps^ ^0, 


(34) 


both in the form of cubic polynomials in z"^. Hence they can be solved analytically, yielding 
bulky explicit formulae for the critical curve and Jacobian saddle points. We proceed to 
identify the boundaries of parameter-space regions with different critical-curve topologies 
using the resultant method described in § 3.2[ applied here to the two cubic polynomials. 


The first resultant yields the polynomial equation 


Pres{w) = (1 — 9fi) -f 6 (1 — 15fi -|- 18fi^) s'^w'^ — 3 (5 J- 3fi) s'^w -|- 8 s® = 0 , (35) 

where w = The second resultant can be factorized to yield three independent condi¬ 

tions. 


8 s® - 3 (5 + 3/x) + 6 (1 - 15/i + 18f?) + 1 - 9fi = 0 

8 s® + 3 (5 + 3/r) + 6 (1 - 15/i + 18f/) - 1 + 9/i = 0 

64 _ 48 (1 - 15/i + 18fi‘^) s® - 3 (5 + 3/i)(l - 9/i) / - (1 - 9fi)^ = 0 


(36) 

(37) 

(38) 
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each defining a curve in parameter space. The left-hand side of equation (36) is equal to 


PresiX)^ hence it describes the splitting of the critical curve at 0 = 0. Similarly, the left-hand 


side of equation (37) is equal to Pres{—Xi hence it describes the splitting of the critical curve 
at 0 = 7r/2. These parameter-space curves thus are genuine boundaries between regions 


with different critical-curve topology. Due to the form of equation (33) any splitting points 


of the critical curve along the real or imaginary axes correspond to parameter combinations 
lying on one of these two boundaries. 


The third equation (38) describes the splitting of the critical curve at off-axis points 


with any other value of 0. Checking the third curve following Appendix we find that for 
/i < 1/5 it is a spurious solution of the resultant method. Only the /i > 1/5 part of the 
curve defines a genuine parameter-space boundary between different topologies. 


m 


A closer inspection of equations (36) and (37) reveals that both are of third degree 
/ while equation (38) is of third degree in s^. Hence, the scale parameter s of all 
the boundaries can be expressed analytically as a function of p. Instead, we present the 
results graphically in the parameter-space plot in the left panel of Figure]^ The boundaries 


described by equations (36), (37), and (38) are illustrated by the black, orange, and cyan 


curves, respectively. The boundaries reach the limits of the plot at [/r, s] = [0,8“^^^], [0,1], 
[1/9,0], and [1,1]. At the [1/5,5“^/^] triple point the critical curve passes through the 
monkey saddles. 


Critical-curve topologies — the plot of the parameter space in the left panel of Figure 
identifies regions with different critical-curve topology. The right panel and its blown-up 
detail in Figure include further subdivision by total cusp number, as discussed further 
below. For better orientation we mark the six topology regions in the left panel by letters 
A - F from top to bottom and left to right. Starting from the widest regime we proceed 
gradually to lower values of the separation s and mark regions in order of appearance. In 
case several regions appear at the same value of s, we mark them in order from lowest to 
highest central mass p. 


The critical-curve topologies corresponding to all regions are illustrated by the examples 
in Figure]^ where they are shown together with their caustics. Figure [^includes examples 
of critical curves and caustics for transitions along the topology boundaries, marked by the 
letters of adjacent topology regions. The last two examples correspond to the boundary- 
intersection points seen in Figure]^ 

As shown in Figure]^ the critical curve of the LS model consists of one (region B), three 
(regions A, C, E), or five loops (regions D, F). Not all the regions have a unique critical-curve 
topology: regions C-l-E, as well as D-l-F have the same topology. The linear symmetric triple 
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lens thus has four distinct topologies: a single loop (B), three independent loops (A), an 
outer loop with two inner loops (C, E), or an outer loop with four inner loops (D, F). The 
topologies are summarized in the LS row of Table 

We inspect the curves in Figure together with the parameter-space map in the left 
panel of Figure in order to understand the occurrence of the different topologies and tran¬ 
sitions between them. In the wide limit (region A) the three independent loops correspond 
to Einstein rings of the three components. Reducing the separation Erst leads to a simulta¬ 
neous merger of all three loops at the real-axis saddles forming a single loop in region B, as 
shown by transition AB in Figure 

Further development with decreasing s depends on the value of fi. For /i < 0.0753 
the next transition BC leads to region C, splitting off inwards two loops around the inner 
pair of Jacobian maxima along the imaginary axis, followed by transition CF to region F, 
splitting from the outer loop another pair of loops around the outer pair of maxima along 
the imaginary axis. For /i G (0.0753,1/9) and /i G (1/9,1/5) the transition BE to region E 
splits off inwards two loops along the imaginary axis, each of which encloses two maxima. 
For fi G (0.0753,1/9) the next transition EF to region F splits each of the inner loops into 
a vertical pair along the imaginary axis, while for /r G (1/9,1/5) transition DE to region D 
splits the inner loops into horizontal pairs bracing the imaginary axis. Finally, for /r > 1/9 
transition BD to region D splits off simultaneously four loops from four different points along 
the outer loop. 

The three remaining special cases, fi ~ 0.0753, /x = 1/9, and /i = 1/5, are illustrated and 
discussed together with their caustic structure under the “Jacobian contour plots” heading 
further below. 

Caustic structure — the number of cusps on the caustic changes in beak-to-beak, butterfly, 
and swallow-tail metamorphoses. These can be found by studying the intersections of the 
cusp and morph curves, as discussed in § |^ The polynomial form of both curves for the LS 
model is provided in Appendix 

The butterfly and swallow-tail metamorphoses introduce additional boundary curves 
that lead to a hner division of parameter space by total cusp number. The number is color- 
coded in the right panel of Figure and in Figure Topology regions A and D are divided 
into two, B into three, and E into four subregions, with regions C and F undivided. For 
better orientation we mark each subregion by the letter of the region and a subscript number 
assigned similarly within the region from top to bottom and left to right. 

A gallery of caustics corresponding to all 13 subregions can be found in Figure [^together 
with transition caustics BCEF and BDE. For each marked subregion there is a pair of panels 
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showing the full caustic plus a blown-up detail. Already at hrst inspection we can identify 
a range of local features that do not occur in two-point-mass lens caustics, such as self- 
intersecting loops (e.g., A 2 , E 2 ), nested (Ci) or overlapping (D 2 , Fi) loops. The caustic 
structures generated in butterfly metamorphoses can be seen in A 2 , B 2 , Di, Ei, E 2 ; swallow- 
tail-generated structures in B 3 , E 4 . 


Inspecting the caustic changes along parameter-space boundaries, we find there are two 
types of beak-to-beak metamorphoses. The transition from two tangent folds to two facing 
cusps may lead to the formation of an additional loop. This is always the case in two-point- 
mass lenses, and it occurs here in transitions B-A (Bi-Ai, B 2 -A 2 ), B-D (B 2 -D 1 ), B-E (B 2 -E 1 , 
B 1 -E 3 , B 3 -E 4 ) and C-F (Ci-Fi). However, we hud the transition may just as well result in 
the opposite, i.e., the merger of two loops. This can be seen here in transitions C-B (C 1 -B 3 ), 
D-E (D 2 -E 2 ), and F-E (F 1 -E 4 ), always starting from overlapping loops and ending in a self- 
intersecting loop. Hence, in the triple lens the formation of a cusp pair in the beak-to-beak 
metamorphosis may either increase or decrease the number of caustic loops. 

Studying the new boundaries in the right panel of Figure]^ we note they are formed by 
three curves, all ending at [/i, s] = [1/9, 0]. The curve leading from [1, \/2 + l] through regions 
A, B, E, and the curve leading from [1, \/2 — 1] through regions D, E are both associated 
with butterfly metamorphoses. The specihc transitions include A 1 -A 2 , B 1 -B 2 , E 3 -E 1 along 
the hrst, and D 2 -D 1 , E 1 -E 2 along the second curve. The third curve leads from [0,1] and 
shadows the adjacent topology boundary from the right side through regions B and E, as 
seen better in Figure This curve is associated with swallow-tail metamorphoses (B 1 -B 3 , 
E 3 -E 4 ). 


The total cusp number in the LS model varies from 8 (subregion Bi marked gray in 
Figures and 1^ via 12 (Ai, B 2 , Ci, E 3 - red) and 16 (A 2 , B 3 , D 2 , Ei, Fi ~ cyan) to 20 
(Di, E 2 , E 4 - brown). Due to the high symmetry of the model, single metamorphoses never 
occur - always at least two occur simultaneously. As a result, the total cusp number in this 
case always changes in multiples of four, with the largest number of cusps appearing in the 
local transition from Bi (8 cusps) to E 4 (20 cusps) just to the right of boundary intersection 
BCEF. There is more diversity in cusp numbers of individual caustic loops, which have 3, 4, 
6 , 8 , 12 (B 2 ) or 16 (B 3 ) cusps. The occurrence of different loop combinations is summarized 
in column LS of Table [2l 


Jacobian contour plots — following the examples given by Danek fc Heyrovsky| (2015), we 
present in Figure [^Jacobian contour plots for three special cases: /i ~ 0.0753, p, = 1/9, and 
p = 1/5. The panels in the upper row illustrate the critical-curve sequences along vertical 
p = const, cuts in the parameter-space plots. In each panel we mark the s = 0.5 critical 
curve in bold for orientation. We recall that for lower s values the critical curves correspond 
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to the outer contours plus the inner loops around the Jacobian maxima, while for higher s 
values they correspond to the inner contours around the lens positions. Note that the central 
(/i = 1/9) and right (/i = 1/5) columns are the only cases of the LS model with four saddles 
instead of the generic six, as discussed above. 

The panels in the lower row include the corresponding cusp (orange) and morph (green) 
curves, which identify the number and distribution of cusps on the corresponding caustics. In 
all cases the contours closest to the lenses have four intersections with the cusp curve around 
each lens, corresponding to three four-cusped caustic loops in the wide-limit subregion Ai. 
The outermost contours corresponding to the close-limit Einstein ring have four cusp-curve 
intersections too. Hence, the central caustic loop has four cusps as well, whether in subregion 
Fi (left), E 3 (central), or D 2 (right column). 

At ^ ^ 0.0753 (left column) the critical-curve sequence is A-B-F, and the caustic se¬ 
quence is A 1 -B 1 -B 3 -F 1 . The transition from B via the boundary intersection BCEF directly 
to region F splits off simultaneously all four inner loops along the imaginary axis. The tiny 
loops around the inner pair of maxima are barely discernible just interior of the inner pair of 
saddles along the imaginary axis (see also Figure]^. Each of the small contours around the 
maxima has three cusp-curve intersections, implying a set of four three-cusped caustic loops 
in the corresponding region Fi. Inspecting the morph curve, note the somewhat hard-to-see 
four simultaneous B 1 -B 3 swallow-tail metamorphoses. These are indicated by the four simple 
intersections of the small figure-eight parts of the morph curve with the small loops of the 
cusp curve. 


In the fi = 1/9 case (central column), the critical-curve sequence is A-B-E, and the 
caustic sequence is A 1 -B 1 -E 3 . Here there are only two double maxima, so that the critical 
curve stays in region E even in the limit s —)■ 0. For this value of fi the close-limit topology 
of the critical curve has an outer Einstein ring with only two small loops inside. In the lower 
panel we see that either of these loops has four intersections with the cusp curve. This means 
that the corresponding caustic loops have four instead of three cusps in the s —)■ 0 limit in 
subregion E3. This configuration is an example of the “multiple caustics” case studied by 


Bozza (2000b) 


The fi = 1/9 case is generally the only case with neither butterfly nor swallow-tail 
metamorphoses, as seen also from the right panel of Figure The cusp and morph curves 
intersect only at the positions of the lenses, saddles, and higher-order maxima. The higher- 


order maxima are isolated solutions of the morph-curve equation (28). When perturbed to 


lower fi values, these points turn into the horizontal figure-eight loops of the morph curve, 
as seen in the left fi ~ 0.0753 column. Similarly, when perturbed to higher fi values, they 
turn into vertical figure-eight loops that eventually merge with the outer part of the morph 
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curve at the /i = 1/5 monkey saddles. 

At /i = 1/5 (right column) the critical-curve sequence is A-B-D, and the caustic sequence 
is A 1 -B 1 -B 2 -D 2 . Here the transition from B to D via the pair of monkey-saddle points at 
BDE splits off simultaneously all four inner loops around off-axis maxima, with the splitting 
occurring pairwise from the two monkey saddles (see also Figure [^. The small contours 
around the maxima each have three cusp-curve intersections, hence they correspond to four 
three-cusped caustic loops in subregion D 2 . Inspecting the morph curves in the lower panel, 
we hrst notice the characteristic structure at the monkey saddles. In each of them the cusp 
curve arrives from four perpendicular directions, while the morph curve arrives from six 
symmetric directions. The corresponding caustic has two beak-to-beak metamorphoses with 
each involving three beaks. These can be seen in the BDE case triple-loop contact in the 
last panel of Figure Note also the additional tangent intersections of the cusp and morph 
curves on the imaginary axis closer to the origin. These correspond to two simultaneous 
B 1 -B 2 butterfly metamorphoses, as expected from Figure 

Note that yet another LS-model Jacobian plot (for the /i = 1/3 equal-mass case) can be 
found further in the right column of Figure where it corresponds to the LA model with 
position parameter p = 1/2. Its LS-model sequence of caustic subregions is A 1 -B 1 -B 2 -D 1 -D 2 . 

Planetary limits — the LS model has two different planetary-mass limits. In the /i —)■ 0 case 
we have an equal-mass binary with separation 2 s and a planet of mass /i placed exactly 
between its components (at the LI Lagrangian point). In order of decreasing s the caustic 
sequence is A 1 -B 1 -B 3 -C 1 -F 1 . At both limits, the planet merely adds tiny extra loops to the 
binary caustic: the central loop in Ai, and the pair of extra three-cusped loops in Ci and Fi. 
The planet perturbs primarily the s = 1 wide - intermediate binary transition, which occurs 
at the position of the planet. Here, in the single-loop subregions Bi and B 3 the planet adds 
extra cusps to the caustic perpendicularly to the binary axis. 


In the fi ^ 1 case we have a star with two equal-mass planets at opposite sides at a 
separation s. Here the caustic sequence is A 1 -A 2 -B 2 -D 1 -D 2 . In the limiting subregions Ai 
and D 2 the planets perturb the point-like single-lens caustic to form the central four-cusped 
loop (Chang & Refsdal 1984), adding two (Ai) or four (D 2 ) extra tiny loops. However, in 
subregions A 2 , B 2 , Di covering the separation range s E [\/2 — 1, \/2 + l] the planets generate 
extra cusps on the central loop at positions perpendicular to the orientation of the system. 
The planets affect the caustic most prominently when they lie in the vicinity of the s = 1 
Einstein radius, causing all loops to merge and form the single-loop caustic B 2 . 


Close limit — for fi 1/9 (subregions Fi, D 2 ): critical curve = Einstein ring -|- four small 
loops around Jacobian maxima; caustic = central four-cusped loop -|- four three-cusped loops 
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escaping to oo. For /r = 1/9 (E 3 ): critical curve = Einstein ring + two small loops around 
Jacobian maxima; caustic = central four-cusped loop + two four-cusped loops escaping to 
also Bozza 2000b( ). 


00 see 


Wide limit — for any ji (subregion Ai): critical curve = three independent Einstein rings 
with radii {^(1 — /x)/ 2 , ^(1 — /i)/ 2 }; caustic = three four-cusped weak-shear Chang- 

Refsdal loops. 


5.3. LA Model: Linear Asymmetric Configuration 


Model description — in the LA model three equal-mass lens components lie along a line with 
a variable separation of the outer two components, marked 2 s to conform with the previous 
LS model. Here the second parameter is the relative separation p of the central component 
from the left component in units of the outer component separation (see second sketch in 
Figure]^. Placing the origin at the geometric center and the lenses along the real axis, the 
positions and masses of the components are {zi, Z 2 , ^ 3 } = {—p — 1, 2p — 1,2 — p} 2 s/3 and 
{pi, /i 2 , ps} = {1/3,1/3,1/3}, with parameter ranges p G (0,1) and s G (0, 00 ). 


The parametrization is symmetric about p = 1/2, so that the conhgurations, critical 
curves, and caustics for p and 1 —p are mutual mirror images. In the p —?• 0 limit components 
1 and 2 merge, so that the lens is reduced to a binary lens with masses (2/3,1/3}. Similarly, 
in the p —1 limit components 2 and 3 merge. For p = 1/2 the components are spaced 
symmetrically, hence these conhgurations are identical to the /i = 1/3 conhgurations of the 
LS model and with the 6 = ti conhgurations of the TI model discussed further in 


5.5 


Jacobian surface character — for any p: four simple maxima -|- six simple saddles. 


Topology boundaries — due to the components having equal masses, we get the center of 
mass di = 0, and the remaining three moments are C 2 = 4s^(l — p + p^)/3, (J 2 = 2 c 2 / 3 , and 
C 3 = 8 s^(—p — l)(2p — 1)(2 — p)/27. The critical-curve polynomial 


Pcrit{z) = - (2 C2 + z‘^ - 2 c 3 Z^ + clz‘^ + 2 03(02 - -F C3 - C2 e^‘'^/3 = 0 (39) 


and the saddle-point polynomial 

Psaddiz) = Z^ + 02Z‘^ + 7o3Z^ + oj z"^ - C3 C2 Z + C3 - 0^/3 = 0 (40) 


generally cannot be solved analytically. We proceed to identify the boundaries of parameter- 
space regions with diherent critical-curve topologies using the resultant method described in 


3.2 applied directly to polynomials Pcru and Psadd- 
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The first resultant yields the polynomial equation 

Pres{w) = (8 — z/) + 6 (12 — u) C 2 + 216 c\ w'^ — (216 — 198 v + 11 u‘^) cl 

—3 u (45 — 14 u) cl — 18 C 2 w + z/^ c® = 0 , (41) 

where w = and the parameter v = 27p^(l —p)^/(l —p + The equation depends on 
s through the moment C 2 . The second resultant can be factorized to yield three independent 
conditions. The hrst, Pres(l) = 0, describes the splitting of the critical curve at 0 = 0 as a 6th 
degree polynomial in s^. The second, pres(—1) = 0, would describe the splitting of the critical 
curve at 0 = 7r/2. However, this condition has no solution in our parameter space, hence 
such splitting does not occur in the LA model. The third condition describes the splitting 
of the critical curve at any other value of 0. It has the form of a 15th degree polynomial in 
which we do not present here explicitly. An inspection following Appendix reveals an 
entire spurious branch of the corresponding curve in parameter space. 

The Pres(l) = 0 and the 15th degree polynomial boundaries are marked in the parameter- 
space plot in the left panel of Figure]^ by the black and cyan curves, respectively. The curves 
meet the p = 0,1 limits at s = — l)“^/^/2 0.981 and s = (^ — l)^/^/2 0.357, 

respectively. 

Critical-curve topologies — the division of parameter space by critical-curve topology is 
shown in the left panel of Figure In view of the p -f-)- (1 — p) symmetry, we use letters 
to mark the six different regions within the left half of the parameter space (from top to 
bottom and left to right), and mark their symmetric counterparts by the same letters. 

The topologies corresponding to all regions are illustrated by the examples presented in 
Figure [^together with their caustics. The LA model permits a greater range of critical-curve 
topologies than the LS model, with each region having a unique topology. In addition to all 
the topologies seen previously in the LS model (Figure]^, there is a two-loop (region B) and 
a four-loop (region A) topology. In region B one of the two independent loops contains two 
lens components, while in region A there is a simple loop plus a close-binary set of loops (two 
small loops within an outer loop). Neither of these can be achieved in the symmetric LS 
model. Note also that the single-loop topology (region D) occurs only for p G (0.239, 0.761). 
In other cases with decreasing s the loop around the pair of components undergoes the 
intermediate - close binary transition (B-A) before connecting with the third-component 
loop (A-E). All topologies are summarized in the LA row of Table 

Caustic structure — the further subdivision of parameter space according to total cusp 
number is shown in the right panel of Figure]^ The number of cusps on the caustic changes 
primarily in beak-to-beak metamorphoses. Unlike in the case of the LS model, all are of 
the standard type with cusp formation accompanying the splitting of caustic loops. Other 
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metamorphoses occur only in the limited range p G (0.292,0.708), as seen from the extent 
of the additional duck-foot-shaped boundary in Figure The curve corresponds to a pair 
of swallow-tail metamorphoses, with the exception of its two points on the p = 1/2 midline. 
The higher symmetry at these points leads to a pair of butterfly metamorphoses at each. 
The new boundary divides regions D, E, and F into two subregions each, while the other 
regions remain undivided. In comparison with the other models, each of the nine subregions 
has a unique caustic structure, as seen from column LA of Table 

The LA model has a lower degree of symmetry than the LS model, so that beak-to-beak 
metamorphoses along the real axis may occur singly and the total cusp number thus may 
change by 2. The total cusp number varies from 8 to 20 in steps of two, with one exception: 
there is no eighteen-cusped caustic in the LA model. The newly occurring totals are 10 cusps 
(subregion Bi marked blue in Figure]^ and 14 cusps (Ai - green). In comparison with the 
LS model, there are caustic loops with 10 cusps (Bi), but there are no loops with 16 cusps. 


Jacobian contour plots — in Figure 10 we present Jacobian contour plots for three cases: 
p = 0.18, 0.45,1/2. The latter two were used by Danek & Heyrovsky (2015) to illustrate the 
swallow-tail and butterfly metamorphoses, respectively. The p = 1/2 caustic sequence from 
wide to close limit is C 1 -D 1 -D 2 -F 1 -F 2 , with transitions D 1 -D 2 and F 2 -F 1 each involving a pair 
of simultaneous butterfly metamorphoses. The p = 0.45 caustic sequence is C 1 -B 1 -D 1 -D 2 -E 2 - 
E 1 -F 2 , with transitions D 1 -D 2 and E 1 -E 2 each involving a pair of simultaneous swallow-tail 
met amorphoses. 


The p = 0.18 case in the left column illustrates a hierarchical triple system, with com¬ 
ponents 1 and 2 forming a binary and component 3 as a distant companion. The structure 
of the contours and curves near the hrst two components resembles a simple equal-mass 
binary (see Danek & Heyrovsky 2015) with its wide, intermediate, and close regimes. If we 
ignore the local structure around components 1 and 2, on the large scale the contours and 
curves again resemble a binary system with its three regimes, here with a 2:1 mass ratio since 
components 1 & 2 act as a single body. The caustic sequence C 1 -B 1 -A 1 -E 1 -F 2 undergoes no 
swallow-tail, no butterfly, and no single-loop topology. 


Planetary limits — the LA model has no planetary limit. 

Close limit — for any p (subregion F 2 ): critical curve = Einstein ring -|- four small loops 
around Jacobian maxima; caustic = central four-cusped loop -|- four three-cusped loops 
escaping to 00 . 


Wide limit — the wide limit is more interesting, with three different regimes for arbitrarily 
large separation s —)■ 00 . In this limit the model behaves as an equal-mass binary lens plus an 
independent single lens. The three regimes are simply the close (subregion Ai), intermediate 
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(Bi), and wide (Ci) regimes of the binary lens. The wide regime, which eventually dominates, 
leads to three independent single lenses with Einstein radii 

This behavior can be proved by the s —>■ oo asymptotic form of the boundaries in 
Figure in the p < 1/2 part the A-B boundary is given by 2ps = -y/I/S, and the B-C 
boundary by 2ps = y/S/S. Noting that 2ps is the separation between lens components 1 
and 2 (see Figure [^, if we divide the two values by the Einstein radius of the binary (y/ 2/3), 
we get the close - intermediate and intermediate - wide boundaries of the equal-mass binary 
lens. To summarize: critical curve = binary critical curve -|- Einstein ring with radius y/ 1/3 ; 
caustic = binary caustic -|- four-cusped weak-shear Chang-Refsdal loop. 

It is worth noticing that the presence of the distant third body along the binary axis 
has no effect on the critical curve topology or caustic structure of the binary lens, it merely 
adds a distant loop. Similarly, the structure of the plots in Figure near the left boundary 
indicates that a tiny separation of components 1 and 2 along the axis to component 3 has 
no effect on the critical curve topology or caustic structure of the 2:1 mass-ratio binary lens, 
it merely adds two tiny loops. 


5.4. TE Model: Triangular Equilateral Configuration 


Model description — in the TE model the lens components he at the vertices of an equi¬ 
lateral triangle with side length s (see third sketch in Figure [^. We use the fractional 
mass /i of component 1 as the second parameter. Placing the origin at the geometric cen¬ 
ter and aligning the real axis with the median passing through component 1, the positions 
and masses of the components are {zi,Z 2 ,zs} = {— 1, s/a/3 and {pi,/i 2 , ps} = 

{ p , (1 — p)/2, (1 — /i)/2}, with parameter ranges p G (0,1) and s G (0, cxd). 


The p = 0 limit corresponds to an equal-mass binary formed by components 2 and 3, 
while the p = 1 limit corresponds to component 1 as a single lens. The p = 1/3 case is the 
most symmetric triple lens: three equal masses in an equilateral conhguration. It is identical 
with the 6* = 7r/3 case of the TI model discussed further in 


5.5 


Jacobian surface character — for /i / 0.768, 8/9: four simple maxima -|- six simple saddles; 
for p ^ 0.768: two simple maxima -|- one double maximum -|- Eve simple saddles; for 
p = 8/9: four simple maxima -|- four simple saddles -|- one monkey saddle. In addition, note 
that the p = 1/3 equal-mass case is unique due to its global three-fold symmetry. 


Topology boundaries — due to the equilateral conhguration C 2 = 0, the center of mass 








di = (1 —3/i) s/\/l2, C 3 = and (i 2 = —(1 —3/i) s^/ 6 . The critical-curve polynomial 


Pcrit{z) = z^-2{c3 + die^^^)z^-3d2e^^^ z‘^-2c3e^^^ z + Csics-die^^^) = 0 (42) 


and the saddle-point polynomial 

Psadd{z) = -h 3 diz^ + 6 d 2 Z^ + 7 Csz^ + Gcsdiz'^ + 3 c 3 d 2 Z + cl = 0 


(43) 


cannot be solved analytically for a general p. We use the resultant method described in § 3.2 


applied to polynomials Pcru and Psadd to identify the boundaries of regions with different 
critical-curve topologies. 


The hrst resultant yields the polynomial equation 


Pres{w) = (r^-3r^-M)w® 

-3r (3r^-12r^ + 8r2+3r-3) sW - 3 r (23r3-20r2-8r + 7) 
-|-(6r^-l-12r^ —11) -|- 6 r (r —2) s^w’^ -f 3t = 0 , (44) 


where w = and the parameter r = (3/i — l)/2 is equal to zero for an equal-mass 
lens. The second resultant can be factorized to yield four independent conditions. The first 
two, Pres(l) = 0 and pres{—7) = 0 , describe the splitting of the critical curve at 0 = 0 
and (j) = 7 r/ 2 , respectively, in the form of 6 th degree polynomials in s^. The other two 
conditions describe the splitting of the critical curve at other values of 0. Both have the 
form of polynomials in (3rd and 12th degree, respectively), which we do not present here 
explicitly. An inspection following Appendix [^reveals that only a segment of the 3rd degree 
polynomial curve corresponds to topology transitions, the rest is spurious. 


The Pres{7) = 0, Pres{—7) = 0, the 3rd, and 12th degree polynomial boundaries are 
marked in the parameter-space plot in the left panel of Figure [TT] by the black, orange, cyan, 
and green curves, respectively. The boundaries reach the p, = 0 side of the plot at s = ■\/l/2 
and s = 2, the p = I side at s = 1, and the s = 0 side at /i = 1 — (4/3) sin (7r/18) ~ 0.768. 
At the [/i, s] = [8/9, ■\/l/ 2 ] triple point the critical curve passes through a monkey saddle. 


Critical-curve topologies — the division of parameter space by critical-curve topology is 
shown in the left panel of Figure H The right panel and its blown-up details in Figure [T^ 
include the further subdivision according to total cusp number. The four curves identified 
above carve the parameter space in the left panel into ten different topology regions, marked 
by letters from top to bottom and left to right. 


For each region examples of critical curves with their corresponding caustics are pre- 
We see that the critical curve may have anywhere from one to hve 


sented in Figure 13 


loops. In comparison with the previous models, there are two new topologies: a small loop 
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within an enter loop (region D), and three small loops within an enter loop (I, J). Only one 
of the previons models’ topologies does not occur here: the simple loop plus a close-binary 
set of loops (model LA region A). As seen from the TE row of Table this model permits 
altogether seven different critical topologies, since the pairs of regions E-l-H, F-l-G, and I-l-J 
each share the same topology. 

The beak-to-beak metamorphoses along most of the topology boundaries increase the 
number of critical-curve loops, with a decrease occurring only for the D-C, I-E, G-J, and 
F-J transitions. In addition, the TE model has the special case of the D-B transition, in 
which two simultaneous beak-to-beaks convert an outer-plus-inner loop combination to two 
separate loops, thus preserving the total number of loops. 


Caustic structure 


the subdivision of parameter space by number of cusps on the caustic 

As seen in both figures, many 


is shown in the right panel of Figure 11 and in Figure 12 


additional metamorphoses change the cusp number. There are four curves corresponding to 
butterfly and one curve corresponding to swallow-tail metamorphoses. 


Two of the butterfly curves start from [0, ^Jl/2]■. one leading down through G to [1/3, 0], 
returning back through G, J, and F to s = —1 -1- a/3 — ■\/2 — a/3 ~ 0.214 at the right edge; 
the second goes up through B and A to s = 1 -|- a/3 -|- ^/2 + a/3 ~ 4.66 at the right edge. 
The other two butterfly curves rise monotonically from [0.768,0] toward the right side: the 
third through J and F to s = 1 -|- y/3 — \/2 + a/3 ~ 0.800; the fourth through J, E, G, and 
A to s = -1 + V3 + ^ 2 - 73 ^ 1.25. 


Finally, the swallow-tail curve starts from [0, ^Jl/2] and leads up through B with a small 
sharp peak at /x ~ 0.084 to A, where it peaks sharply when it meets the butterfly curve at 
fi = 1/3, dropping down through A, G, E, and J to [0.768,0]. The curve corresponds to a 
complex-conjugate pair of swallow-tail metamorphoses everywhere except at the two sharp 
peaks. At the /i = 1/3 peak three symmetric butterfly metamorphoses occur simultaneously, 
while at the ja ~ 0.084 peak there is a complex-conjugate pair of butterfly metamorphoses. 


The additional set of curves in the right panel of Figure 11 (with details better seen in 
Figure divides regions G, E, F, and G into three subregions each, regions A and B into 
five subregions, and region J into eight subregions. Regions D, H, and I remain undivided. 
The parameter space of the TE model is thus divided into 33 subregions, corresponding to 
24 different loop combinations. These are sorted by cusp number and listed in column TE 
of Table The total cusp number varies from 10 to 20 in steps of two; i.e., there is no 
eight-cusped caustic. Gaustics with 18 cusps that appear newly in this model occur in seven 
different subregions (see Table [^. Individual caustic loops have 3, 4, 6 , 7, 8 , 9, 10, 12, or 
14 cusps. In comparison with the previous models there is no sixteen-cusped loop (B 3 in LS 
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model), while seven-cusped (Ii shown in Fignrej^ and nine-cusped loops (Di in Fignre[I^ 
are new here. 


Jacobian contour plots — Fignre 14 inclndes contonr plots for three special cases: /i = 1/3, 
pi ~ 0.768, and p. = 8 / 9 . In the symmetric p = 1/3 case the canstic seqnence from 
wide to close limit is A1-A4-D1-G2, with the A1-A4 three simnltaneons bntterfly meta¬ 
morphoses at s = 2^/^ 3^/^ ~ 2.75. The D1-A4 three simnltaneons beak-to-beaks occur 

at s = [(5\/5 -|- ll)/2]^/® 1.49, and the D1-G2 three simultaneous beak-to-beaks occur 

at s = [(5a/5 — ll)/2]^/® 0.670. There is no single-loop topology, but with the slight¬ 

est perturbation to a higher p the sequence would visit the G2 single-loop region with a 
fourteen-cusped caustic. 


The double-maximum case in the central p ~ 0.768 column has a caustic sequence Ai- 
A2-G1-E2-J3-J7 with two butterfly metamorphoses along the real axis. The plot requires a 
larger scale due to the close-limit J7-J3 butterfly point near the right edge of the inset. 

The same holds also for the p = 8/9 monkey-saddle case in the right column. The 
caustic sequence A1-A2-G1-G3-E1-F2-F3 starts out similarly, but differs once the components 
get closer. The three butterfly points all lie on the real axis. 


Planetary limits — the TE model has two different planetary-mass limits, similar to those 
of the LS model. The /i —)■ 0 limit describes an equal-mass binary separated by s with 
a planet of mass p offset perpendicularly at the vertex of an equilateral triangle (at the 
L4 Lagrangian point for a face-on coplanar orbilj^. In order of decreasing s, the caustic 
undergoes a sequence of seven metamorphoses: A1-B1-B5-B3-D1-II1-G2-G1. The two widest 
regimes, Ai and Bi, correspond to a wide and intermediate binary caustic, respectively, each 
with an additional four-cusped loop due to the planet. Similarly, Gi corresponds to a close 
binary caustic with two additional three-cusped loops due to the planet. 

All six metamorphoses between Bi and Gi occur close to the s = 1/a/2 intermediate 
- close binary transition, since the planet lies exactly at the location at which one three- 
cusped loop would split off in the binary limit. The sequence begins with a butterfly on the 
planetary loop leading to B5, followed by a pair of swallow tails on the binary-caustic loop 
to B3. In the peculiar transition to Di both loops connect in a pair of reverse beak-to-beaks, 
which simultaneously split off a three-cusped loop. A second three-cusped loop is formed 
along the real axis in a beak-to-beak to Hi, and a third and fourth three-cusped loop appear 
simultaneously in a pair of beak-to-beaks to G2. Finally, the main caustic loop loses two 
cusps in a reverse butterfly to Gi. 


^Note that in this case L4 is unstable, since y .2 = Us- 
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The jU —>■ 1 limit describes a single lens with two equal-mass planets at the same distance 
s, an angle of vr/S apart. Even in this case the caustic undergoes seven metamorphoses with 
decreasing s: A1-A2-A5-C3-E1-F1-F2-F3. In the Ai wide limit, the caustic consists of three 
four-cusped weak-shear Chang-Refsdal loops ( Chang fc Refsdai||1984 ). In the F 3 close limit, 
the caustic has one four-cusped weak-shear Chang-Refsdal loop for the star, plus two pairs 
of three-cusped strong-shear Chang-Refsdal loops for the planets. 


The butterfly metamorphosis from Ai to A 2 , found on the caustic loop of the star at 
the cusp facing the planets, occurs already at a large distance of the planets - approaching 
s ~ 4.66 for p —)■ 1. The next butterfly to A 5 near s ~ 1.25 occurs on the caustic loop of the 
star at the tip of the hrst butterfly feature. The three following beak-to-beak transitions all 
occur near s = 1, the Einstein radius of the single lens. First, the two planetary four-cusped 
loops connect with the star’s loop leading to C 3 , one pair of three-cusped planetary loops 
then disconnects to Ei, and another pair disconnects to Fi. The last two metamorphoses 
occur on the star’s loop: reverse butterflies near s ~ 0.800 to F 2 , and s ~ 0.214 to F 3 . We 
point out that even though the planetary loops of the caustic connect and disconnect when 
the planets he close to the star’s Einstein ring, the star’s loop of the caustic is affected by 
having extra cusps for a much wider range of planetary distances, s G (0.214, 4.66). 

Close limit — for /r 7 ^ 1/3,0.768 (subregions F3, Gi, G3): critical curve = Einstein ring 
-|- four small loops around Jacobian maxima; caustic = central four-cusped loop -|- four 3- 
cusped loops escaping to 00. For fi = 1/3 (G2): critical curve = Einstein ring -|- four small 
loops around Jacobian maxima; caustic = central six-cusped loop -|- three three-cusped loops 
escaping to 00 -|- one three-cusped loop staying at center. For fi = 0.768 (J7): critical curve 
= Einstein ring -|- three small loops around Jacobian maxima; caustic = central four-cusped 
loop -|- one four-cusped and two three-cusped loops escaping to 00. 


The fi = 1/3 equal-mass case exhibits a new close-limit behavior: the central caustic 
loop has six instead of four cusps, and one of the four three-cusped loops stays at the center 
instead of escaping to 00 (see Figure 13). The central caustic loop is self-intersecting with 
a three-fold symmetry, so that a pair of cusps points in the direction of each of the lens 
components. In the s —)■ 0 limit each pair rapidly converges to a point, so that the shape of 
the shrinking caustic loop approaches a simple symmetric three-cusped loop, which is traced 
twice as one follows the corresponding Einstein ring. The nature of the central loop can 
also be seen from the contour plot in the left column of Figure 14 The six branches of the 
cusp curve extending symmetrically to the outer edges have six intersections with the outer 
contours, which represent the perturbed Einstein-ring loop in the close limit. 


Wide limit — for any fi (subregion Ai): critical curve = three independent Einstein rings 
with radii {y/jl, y/(l — /i)/2, ■y/(l — fi)/2}-, caustic = three four-cusped weak-shear Ghang- 
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Refsdal loops. 


5.5. TI Model: Triangular Isosceles Configuration 


Model description — the final, topologically richest model has three equal-mass lens com¬ 
ponents lying at the vertices of an isosceles triangle, with legs of length s spanning a vertex 
angle 6 (see fourth sketch in Figure[^. If we place the origin at the geometric center, align the 
axis of symmetry with the real axis, and set component 1 at the vertex angle, the positions 
and masses of the components are {zi,Z 2 ,Z 3 } = {—-|- 2 2 — 

g-ie/ 2 }_g /3 {jui, jU 2 , jUsj = {1/3,1/3,1/3}, with parameters 6 G (0, vr] and s G (0, cx)). 

In the 9 = 0 limit components 2 and 3 coincide, reducing the system to a binary lens 
with masses (1/3, 2/3} and separation s. We recall that in the LA model both p = 0 and 
p = 1 edges correspond to a binary lens with the same mass ratio but with separation 2 s. 
The 6 = 71 case corresponds to the collinear limit with component 1 at the center. Hence, 
the corresponding conhgurations are identical to the /i = 1/3 case of the LS model, as well 
as the p = 1/2 case of the LA model. Similarly, conhgurations with vertex angle 6 = tt/S 
are identical to the /i = 1/3 conhgurations of the TE model. 


Jacobian surface character — for 6 ^ 0.124 vr, 0.708 vr, 0.847 vr: four simple maxima -|- six 
simple saddles; for 9 k. 0.708 tt : two simple maxima -|- one double maximum -|- hve simple 
saddles; for 9 th 0.124 tt and 9 ~ 0.847 tt : four simple maxima -|- four simple saddles -|- one 
monkey saddle. In addition, the 6 * = 7 r /3 equilateral case has global three-fold symmetry, as 
discussed in 


5.4 


Topology boundaries — since the masses are equal, the center of mass di = 0, C 2 = [1 — 
4 sin^ (6^/2)] s^/3, d 2 = 2 c 2 / 3 , and C 3 = —2 cos (6'/2)[l -|-8 sin^ (0/2)] s^/27. The critical- 
curve polynomial 


Pcrit{z) = z^ — {2c2 + z'^ — 2 C3 z^ c\ z'^ + 2 €3(02 — z + c^ — C2 e^“^/3 = 0 (45) 


and the saddle-point polynomial 

VmddiA = z‘+ C2Z‘+ 7 C3Z^ + - C2C3 Z - C ^/3 + c| = 0 


(46) 


cannot be solved analytically for a general 9. We use the resultant method described in § 3.2 


applied to polynomials pcrit and Psadd to identify the boundaries of regions in parameter 
space with diherent critical-curve topologies. 
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The first resultant yields the polynomial equation 


Pres{w) = [8 (1 “ + 27k] tn® + 6 (1 — k) [4 (1 — k)^ + 9 k] 

+24 (1 - Kfs^w^ - [8 (1 - k)® + 198 K (1 - Kf + 297 n^] 

+9 K (1 - k) [5 (1 - Kf + 42 k] sW - 54 - nfs^^w - 27 = 0, (47) 


where w = and the parameter k = 4 sin^ (^/2) is equal to unity for an equilateral 
conhguration. The second resultant can be factorized to yield three independent conditions. 
The hrst two, Pres(l) = 0 and pres(—1) = 0, describe the splitting of the critical curve at 
0 = 0 and 0 = 7 r/ 2 , respectively, in the form of 6 th degree polynomials in s^. The third 
condition describes the splitting of the critical curve at other values of 0. It has the form of a 
15th degree polynomial in s^, which we do not present here explicitly. An inspection following 
Appendix [^reveals that two segments of the corresponding curve have to be omitted, since 
they yield spurious solutions with no change in critical-curve topology. 

The Pres(l) = 0, Pres(~l) = 0, and the 15th degree polynomial boundaries are marked in 
the parameter-space plot in the left panel of Figure[^by the black, orange, and cyan curves, 
respectively. The boundaries reach the 6 ^ = 0 side of the plot at s = (-^ — 1)^/^ 0.714 

and s = (-0^ — 1)“^/^ 1.96, the 9 = ti side at s 0.546 and s k, 1.68, and the s = 0 side 

aX 9 K. 0.708 TT. At either of the [0,5 ] [0.124 tt, 1.76] and [0.847 vr, 0.478] triple points the 

critical curve passes through a monkey saddle. 


Critical-curve topologies 


shown in the left panel of Figure 15 


the division of parameter space by critical-curve topology is 
The cusp-number map in the right panel and its 


blown-up details in Figure include further subdivision by cusp number. The three curves 
obtained from the resultant analysis carve the parameter space into twelve different regions, 
marked by letters A - L in the left panel, from top to bottom and left to right. 


For each region, examples of critical curves with their corresponding caustics are shown 
in Figure]^ The critical curve may have anywhere from one to hve loops. The TI model has 
all the topologies found in the three previous models, plus one new topology: two separate 
loops, one of which contains a small third loop (region B). As seen from the TI row of 
Table this model permits altogether nine different critical topologies, since the pairs of 
regions E+J, I+L, and H+K each share the same topology. 


The beak-to-beak metamorphoses along most of the topology boundaries increase the 
number of critical-curve loops, with a decrease occurring only for the F-G, I-J, H-L, and K-L 
transitions. Similarly to the D-B transition of the TE model, the TI model has the special 
F-C transition, in which two simultaneous beak-to-beaks convert an outer-plus-inner loop 
combination to two separate loops, preserving the total number of loops. 
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and in Figure [ 16 ] 

include subdivision due to additional caustic metamorphoses. Four of the hve additional 
curves dehne conditions for the butterfly and one for the swallow-tail metamorphoses. Two 
of the butterfly curves asymptotically approach 0 —)■ 0 as s —)■ oo: one leads down through 
A, E, and H to [0, s] = [ 7 r/ 3 , 0 ], rising back through H, L, and K to s ~ 0.502 at the right 
edge; the other leads down through A, B, and C, rising up through C and D to s = oo at 
6 = 7 r/ 2 , dropping down through D and G to s ~ 1.20 at the right edge. The second pair of 
butterfly curves rises from [0.708 7 r, 0 ] to the right edge: one through L and K to [tt, 0.502]; 
the other through L, J, and G to [vr, 1.20]. All the butterfly metamorphoses dehned by these 
four curves occur along the symmetry axis of the lens. 

The swallow-tail curve starts asymptotically at s —)■ cx) as 6^ —)■ 0, dropping down through 
A, B, and G, rising up through G and D to a sharp peak at 6 * ~ 0.220 vr, bouncing back to 
another sharp peak reaching the second butterfly curve at 6 = 7 r/ 3 , dropping hnally through 
D, G, J, and L to [0.708 vr, 0]. Just as in the TE model, the curve corresponds to a complex- 
conjugate pair of swallow-tail metamorphoses everywhere except at the two sharp peaks. At 
the 6 = 71 /3 peak, three symmetric butterfly metamorphoses occur simultaneously, while at 
the 6 ~ 0.220 vr peak there is a complex-conjugate pair of butterfly metamorphoses. These 
are the only two cases within the TI model with butterflies occurring off the symmetry axis. 

The additional set of curves divides region E into two subregions, regions B, G, H, J, 
and K into three subregions each, regions A and G into four, region D into six, and region 
L into eight subregions, with regions F and I remaining undivided. The parameter space of 
the TI model is thus divided into 41 subregions, corresponding to 28 different caustic-loop 
combinations. All are sorted by cusp number and listed in column TI of Table The total 
cusp number varies from 8 to 20 in steps of two. Individual caustic loops have 3, 4, 5, 6 , 7, 
8 , 9, 10, 12 or 14 cusps. In comparison with the previous models there is no sixteen-cusped 
loop (B 3 in LS model), while hve-cusped loops appear newly here (in subregion Bi). 

Jacobian contour plots — Figure [T^ includes contour plots for three special vertical transects 
through the parameter space. Vertex angle 6 ^ 0.220 tt in the left column corresponds to 
the left sharp peak of the swallow-tail curve in Figure The caustic sequence from wide 
to close limit is D1-D4-G2-G3-F1-E2-H2-H1, with no single-loop topology and no swallow-tail 
metamorphosis. The inset includes the II1-II2 butterfly metamorphosis close to its right 
edge. Note the special F-G pair of beak-to-beaks at the Re[z] ~ —0.14 saddle points, in 
which two nested loops reconnect to form two separate loops. 

The D 1 -D 4 simultaneous complex-conjugate pair of butterflies is indicated by the morph 
curve passing through self-intersection points of the cusp curve near Re[z] ~ 0.15. These 
self-intersections disappear for any perturbation of 6 . For lower values the cusp curve at the 


Caustic structure — the cusp-number maps in the right panel of Figure 15 
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upper point splits along a SW-NE axis, for higher values along a NW-SE axis. In either 
case only one of the two separated branches intersects the morph curve, corresponding to a 
complex-conjugate pair of swallow tails. 

The central 6 = 7 r /2 column has the caustic sequence D 2 -G 2 -G 3 -F 1 -I 1 -H 2 -H 1 . The G 2 - 
G 3 pair of swallow tails occurs in the image plane on the central small figure-eight loop of 
the morph curve, and the H 1 -H 2 butterfly can be seen at the left edge of the inset. The 
peculiar feature occurring at this angle is the six-cusped loop around the left component 
in the wide limit (D 2 ). In the double-maximum 6 ~ 0.708tt case in the right column the 
caustic sequence D 3 -D 2 -G 2 -J 2 -L 4 -L 7 ends with a four-cusped caustic loop corresponding to 
the double-maximum critical-curve loop. 


Two further transects of the TI model can be found in contour plots of other models. 
The symmetric 6 = tt/ 3 case appears as the /i = 1/3 left column of the TE-model Figure [l4| 
(TI-model caustic sequence: D 1 -D 6 -F 1 -H 2 ). The linear 9 = tt case appears rotated by 7r/2 
as the p = 0.5 plot of the LA model in Figure 10 (Tl-model: D 3 -G 1 -G 4 -K 1 -K 3 ). 

Planetary limits — the TI model has no planetary limit. 


Close limit — for 6 ^ 7 ^ 7r/3,0.708 tt (subregions Hi, H3, K3): critical curve = Einstein ring -|- 
four small loops around Jacobian maxima; caustic = central four-cusped loop -|- four three- 
cusped loops escaping to 00. For 9 = tt/ 3 (H 2 ): critical curve = Einstein ring -|- four small 
loops around Jacobian maxima; caustic = central six-cusped loop -|- three three-cusped loops 
escaping to 00 -|- one three-cusped loop staying at center. For 9 = 0.708 tt (L7): critical curve 
= Einstein ring -|- three small loops around Jacobian maxima; caustic = central four-cusped 
loop -|- one four-cusped and two three-cusped loops escaping to 00 . 


Wide limit — it is the wide limit that is particularly interesting here. For any fixed 9 the 
s —)■ cxD limit corresponds to three single lenses with Einstein radii a/ 1/3. However, for 
9 ~ the wide limit describes an equal-mass binary formed by components 2 and 3 with a 
distant object offset perpendicularly from the binary axis. This is analogous to the p ~ 
wide limit in the LA model with a distant object lying along the binary axis (see § 5.3). The 


TI model has a total of nine different caustic regimes for arbitrarily large leg length s —>■ 00 . 
In order of increasing vertex angle 9 these are: Ai, A 2 , A 3 , A 4 , Bi, Gi, Di, D 2 , and D 3 . The 
regimes can be grouped by the character of the binary lens formed by components 2 and 
3: close (topology region A), intermediate (G), and wide (D). Topology B and the caustic 
subregion structure are the results of perturbation of the binary by the distant companion 
- as well as perturbation of the companion by the binary. 

The additional topology B (see Figure [ItI) occurs at the close - intermediate boundary. 
The influence of the companion causes the opposite inner loop to connect with the outer 
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loop before the adjacent one, rather than both connecting simnltaneonsly. In the asymp¬ 
totic s —>■ oo regime B becomes negligible, as can be seen from the topology bonndaries. 
The separation of binary components 2 and 3 (see Fignre along the A-B bonndary is 
asymptotically 2s sin( 6 */ 2 ) ~ (1 - 1 - s “^/6 — s“^/6)/\/3, while along the B-C bonndary we 
get 2 s sin (6*/2) ~ (1 -|-s“^/ 6 -|- s“^/6)/a/3. The two expressions differ only at the order. 
For completeness, the asymptotic C-D bonndary yields 2 s sin (0/2) ~ -^8/3. Normalizing 
the asymptotic separations by the ^2/3 Einstein radins of the binary, we recover the stan¬ 
dard close - intermediate A-(B)-C and intermediate - wide C-D bonndaries of an eqnal-mass 
binary. 

The wide-limit canstic in regions A throngh C always has a fonr-cnsped weak-shear 
Chang-Refsdal loop for the distant companion. The binary part of the canstic is dominated 
by the close-binary strnctnre in Ai: three loops with 4-I-3-I-3 cnsps. Increasing the separation 
canses the main loop to nndergo two snbseqnent bntterfly metamorphoses to A 2 (6-I-3-I-3) 
and A 3 (8-I-3-I-3), followed by a pair of reverse swallow tails to A 4 (4-I-3-I-3). One three-cnsped 
loop merges with the main loop to Bi (5-1-3), followed by the other to reach the intermediate 
regime Ci (6), which then splits into two in the wide-binary regime Di (4-1-4). The strnctnre 
of the binary part does not change fnrther in the D region. 

At larger vertex angles in region D, it is the companion loop of the canstic that nndergoes 
an interesting transition. The bntterfly entrance to D 2 tnrns the loop into a six-cnsped self- 
intersecting cnrve. This strnctnre disappears in a reverse bntterfly when exiting to D3. While 
for any other angle an increase in s wonld eventnally retnrn the six-cnsped loop to the regnlar 
fonr-cnsped regime, for 9 = 7 r /2 the pecnliar shape persists for s —)■ 00 . In this limit each 
of the three pairs of cnsps converge to a point, so that the shape approaches a three-cnsped 
loop traced twice. The canstic loop is similar to the six-cnsped loop in the close limit of the 
eqnal-mass TE model. 

We stress that in this regime all three components lie at arbitrarily large separations. 
Still, the conhgnration alters the canstic loop of the component at the vertex fnndamentally. 
This effect is cansed by the gravitational fields of the other two components, which combine 
so that they exactly cancel the lensing shear term. Hence, the Chang-Refsdal limit is not 
valid in this special case. Even thongh this occnrs in the TI model for s —)■ cx) only for 
6 = 7r/2, for any fixed not-too-large valne of s there is an interval of angles, for which the 
simple Change-Refsdal shape of the canstic loop cannot be assnmed. 
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5.6. Overview of critical-curve topologies and caustic structures 

In Table [T] we summarize the occurrence of the nine different critical-curve topologies 
in the parameter spaces of the four studied triple-lens models. The topologies are listed hrst 
in order of increasing total number of loops, then by decreasing number of separate outer 
loops. 

The total loop number ranges from one to five, with single 1-loop and 5-loop topologies, 
two different 2-loop and 4-loop topologies, and three different 3-loop topologies. The number 
of outer separate loops ranges from one to three, and the number of inner loops ranges from 
zero to four. Four of the topologies occur in all studied models: single loop, three separate 
loops (the generic wide limit), an outer loop with two inner loops, and an outer loop with 
four inner loops (the typical close limit). These four are the only topologies found in the LS 
model. The LA model has six, the TE model seven, and the TI model has all nine topologies. 

The topologies found in the analyzed models are only examples of the possible topologies 
of triple-lens critical curves, so the list in Table is by no means exhaustive. All of the 
studied models are symmetric by their geometry and mass combination; all have at least 
one reflection symmetry. Breaking the symmetries of the models opens up the possibility 
of additional critical-curve topologies. However, finding the overall number of topologies of 
the triple lens is beyond the scope of this work. Such a survey would require a detailed 
investigation of the geometry of the Jacobian surface and the positions and types of its 
maxima, saddles, and poles. The correspondence with Jacobian contours may prove to be 
another useful tool, which may identify the possible sequences of topology changes when 
scaling the lens from the wide to the close limit. 

Finally, we point out that the topologies discussed here are those that occur in two- 
dimensional areas of the parameter space of their models. One may proceed further and 
classify also the critical-curve topologies of all transitions between the different marked 
topology regions, such as those shown in Figure While this would be a straightforward 
extension, we do not include it here since the present work is extended enough already. 

Table presents an overview of all 32 different caustic structures in the parameter 
spaces of the four studied triple-lens models. Here we identify the caustic structure purely 
by the combination of cusp numbers on the individual loops. In the notation used in the 
first column of the table, for example 16/10-1-3-1-3 identifies a caustic with a total 16 cusps 
on three loops: one with 10, and two with 3 cusps each. The caustics are ordered in the 
table first by increasing total cusp number, then in decreasing order by the number of cusps 
on individual loops. 

The total cusp number runs from 8 to 20 in increments of 2. The eight-cusped caustics 
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have only a single-loop structure, the ten-cusped caustics have two, twenty-cusped caus¬ 
tics three, eighteen-cusped caustics hve, twelve- and fourteen-cusped caustics six each, and 
sixteen-cusped caustics nine different structures. Individual loops have 3 - 10, 12, 14, or 16 
cusps. Least common are the five-cusped loop (the 12/5-1-4-1-3 caustic in the TI/Bi subre¬ 
gion), and the sixteen-cusped loop (the 16/16 caustic in the LS/B 3 subregion). There are only 
four caustic structures that occur in all four studied models: 12/12, 12/4-1-4-1-4 (includes the 
generic wide limit), 16/4-|-3-|-3-|-3-|-3 (includes the generic close limit), and 20/8-|-3-|-3-|-3-|-3. 
The LS model has 10, the LA model has only 9, the TE model has 24, and the TI model 
has 28 different caustic structures, lacking only four of all the found possibilities. 

The list in Table is based purely on the cusp numbers on loops of the caustic. If we 
took into account the combination of critical-curve topology and caustic structure, the list 
would expand to 37 different cases. The reason for this is that a given caustic structure may 
occur for different mutual positions of critical-curve loops. For example, two four-cusped 
loops of the caustic may correspond to two separate loops of the critical curve, or to an 
outer and inner loop of the critical curve. An inspection of the topologies in Table shows 
that this is the case for the 12/4-1-4-1-4 caustic in the LS model from the A versus E regions; 
the 16/6-1-6-1-4 caustic in the TI model D region versus LS model E region; the 14/4-|-4-|-3-|-3 
caustic in the LA and TI model A regions versus the TE model J region and TI model L 
region; the 16/6-|-4-|-3-|-3 caustic in the TI model A region versus the TE model J region 
and TI model L region; and the 18/8-|-4-|-3-|-3 caustic in the TI model A region versus the 
TE model J region and TI model L region. 

The list in Table does not include the caustics of lens conhgurations at metamorphosis 
points. The classification along subregion boundaries in parameter space would require 
additional counting of beak-to-beak, swallow-tail, and butterfly points. Even here we note 
that the caustic structures of the analyzed models are only examples of triple-lens caustics, 
and the list in Table is not exhaustive. Analyzing less symmetric systems reveals that 
the frequency of metamorphosis types is different from most of the analyzed models, with 
swallow tails occurring more frequently than butterflies (as seen here in the LA model). 


6. Summary 


We presented an initial exploration of the character of critical curves and caustics of 
triple gravitational lenses. The topic is of particular interest for the analysis of observed 
planetary microlensing events. Already four published events involved triple-lens systems: 


in two cases the lens was a star with two planets (( 

daudi et ah 

2008 

Han et ah 

2013 

); in 

two cases a binary lens with a planet ( 

Gould et ah 

2014 

Poleski et ah 

2014 

). Moreover, 
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even triple systems formed by a star with a planet with a moon might be detectable by 
microlensing (e.g., Liebig & Wambsganss 2010|. 


Our approach is based on analyzing the properties of the lens in the image plane, with 
the main relevant results from Danek & Heyrovsky (2015) introduced in §[^ We describe in 
detail in § [^ and § several analytical and numerical methods for parameter-space mapping 
of the topology of the critical curve and the number of cusps of the caustic. The methods 
are presented in a form applicable to any n-point-mass lens. 


In § |5.1 we present the main relevant polynomial equations for triple lenses in a compact 
form in terms of moments of the mass distribution. In § 5^ - §111 we apply the described 
methods to the analysis of four simple two-parameter models of triple lenses. For each of 
them in turn we discuss separately the Jacobian surface character, topology boundaries, 
critical-curve topologies, caustic structure, Jacobian contour plots, planetary limits, close 
limit, and wide limit. The combined results explain the properties of lenses on four 2D cuts 
through the full 5D parameter space of the general triple lens. Since each of the studied 
models intersects at least one other along a ID set of parameter combinations, the presented 
results may serve as a reference framework for further studies. 

The results include the description of a range of triple-lens features, such as the oc¬ 
currence of swallow-tail and butterfly metamorphoses, Jacobian monkey saddles, or double 
maxima. We demonstrated that the appearance of two new cusps in the beak-to-beak meta¬ 
morphosis is not always accompanied by the splitting of critical-curve loops. It may lead to 
the opposite, the merger of loops, as discussed in the caustic structure paragraphs of § |5.2 


Finally, two simultaneous beak-to-beaks may lead to a reconnection preserving the number 
of loops. Within the studied models this occurs along the D-B transition in the TE model, 
and along the F-C transition in the TI model. 


The planetary limits are particularly noteworthy in the equilateral TE model. The 
/i —)■ 0 limit shows that a planet placed near the close - intermediate beak-to-beak point of 
a binary lens may turn the single transition into a sequence of six caustic metamorphoses. 
The /i —)■ 1 limit demonstrates that a pair of planets may form additional cusps on the 
Chang-Refsdal four-cusped caustic of the star at a very broad range of separations, in this 
case s G (0.214, 4.66) Einstein radii. 


Unusual behavior in the close limit occurs in the presence of double maxima (LS, TE, 
and TI models), and in the equal-mass equilateral lens (TE and TI models). Any double 
maximum causes the lens to have a four-cusped escaping caustic loop instead of a pair of 
three-cusped loops (Bozza 2000b). The symmetry of the equal-mass equilateral lens causes 
it to have a self-intersecting six-cusped (asymptotically three-cusped) primary caustic loop. 
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instead of the generic fonr-cnsped loop. In addition, one of the weak three-cnsped canstics 
stays at the origin instead of escaping away. 


Very interesting wide-limit behavior was fonnd in the TI model with vertex angle 
9 = TT12. Here one of the canstic loops retains a self-intersecting six-cnsped shape approach¬ 
ing a three-cnsped shape for arbitrarily large separations of all components. Snch a lens 
violates the generic Chang-Refsdal limit (Chang & Refsdal 1984). In this case the gravita¬ 
tional helds of the two other components exactly combine to cancel the lensing shear term. 
The shape is thns dominated by higher-order terms, leading to its different geometry. 


For smaller angles, the wide limit of the TI model demonstrates the sensitivity of the 
binary canstic strnctnre of the two close components to a distant companion. Here the 
presence of the third perpendicnlarly offset body pertnrbs primarily the canstic strnctnre 
near the close - intermediate transition. As a resnlt, the canstic nndergoes a seqnence of 
five metamorphoses instead of a single one. 


We compiled the resnlts fonnd for all fonr models in the tables presented in § 5.6 The 


critical cnrves have 9 different topologies, and the canstics have 32 different strnctnres, when 
identified by the combination of nnmbers of cnsps on individnal loops. A joint classification 
of critical-cnrve topologies and canstic strnctnres increases the nnmber farther to 37 different 
sitnations. The lists are jnst a sample of the fnll range provided by the general triple lens. 


In Table we provided for each model the parameters of examples of all different canstic 
strnctnres. These can be nsed to guide more detailed studies of critical-curve and caustic 
metamorphoses, to generate examples for testing light-curve fitting codes and algorithms, or 
as a starting point for exploring the behavior of less symmetric triple lenses. 


We thank the anonymous referee for suggestions that helped improve the manuscript. 
Work on this project was supported by Czech Science Foundation grant GACR P209-10-1318 
and Charles University grant SVV-260089. 


A. Computing the resultant of two polynomials 

A polynomial f{x) of degree d>l with roots .^i,... .^^ can be written as 

d d 

/(^)= 

j=0 j=l 
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where 7 ^ 0. Similarly, for another polynomial g{x) of degree e > 1 with roots rji,.. .rjf, we 
have 

e e 

j=0 j=l 

with be 7 ^ 0. The resnltant of the two polynomials is a fnnction of their coefficients that 
is proportional to the prodnct of the differences of all combinations of individnal roots, 


specifically (e.g., Escofier 2001), 


Res^(/, g) = a^dKWW (0 -Vk)- 

j=i k=i 


(A3) 


Using eqnations (Al) and (A2), we can also write the resnltant in the form (e.g., Stnrmfels 

d e 

Res^(/, c/) = alW g{ij) = JJ /(r/fc). 


2002 ) 


(A4) 


i=i 


k=l 


It follows from eqnation (|A3|) that the condition for the two polynomials to have at least one 

(A5) 


common root is eqnivalent to the condition 

Res^(/,5') =0, 

an eqnation in terms of the coefficients of the polynomials. 


A common way to compnte the resnltant is based on the Sylvester matrix S'(/, g). This 
(d + e) X (d + e) matrix is formed by arranging the coefficients of f{x) from highest to lowest 
order, filling the seqnence repeatedly in the first e rows gradnally staggered to the right, 
followed by d rows of coefficients of g{x) staggered in a similar manner, filling the rest with 


zeros (e.g.. Fetters et al. 2001). Explicitly, 

^ Od-l 
0 ad 


S{f,g) = 


ao 0 

Ol Qq 


0 0 
be—1 
0 be 


bo 0 
bi bo 


0 \ 
0 


ao 

0 

0 


(A 6 ) 


y 0 0 ■ • • bo / 

The resnltant is then obtained simply by compnting the determinant of the Sylvester matrix, 

Res^(f,g) = detS(f,g). (A7) 
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In some texts (e.g., Escofier 2001) the Sylvester matrix is defined as the transpose of the 


matrix in equation (A6), i.e., with coefficients in staggered columns rather than rows. How¬ 


ever, this does not change the determinant, so that equation (A7) remains valid. In other 


texts the form of the matrix in equation (A6) remains the same, but the coefficients are 


filled in reverse, from lowest to highest order (e.g., Erdl & Schneider 1993 Sturmfels 2002). 
The determinant of such a matrix is not always equal to the resultant defined by equa¬ 
tion (A3), it differs by a factor (—1)'^®. However, this difference plays no role when seeking 


the null-resultant condition given by equation (A5), which is equivalent to det S{f,g) = 0. 


An alternative way to compute the resultant utilizes the smaller n x n Bezout matrix 
B{f,g), where n = Ma.x{d,e). The matrix is obtained by anti-symmetrizing the direct 
product of the polynomials, dividing the result by the difference of their variables, and 
ordering the terms according to the powers of the variables. The elements of the matrix are 


thus defined by (e.g., Sturmfels 2002) 


n—1 

Bjk{f,g)x^y 

j,k=0 


k _ fix)g{y)-f{y)g{x) 
x-y 


(AS) 


If we assume that the degree of the first polynomial d > e, substituting the expressions from 


equations (Al) and (A2 ) in equation (A8) then yields Bjk{f,g) in terms of the coefficients: 


d-l 


d-l 


Min(d, ji+Ai+l) 

a) X^y^ = y'" Y bj+k-i+i 

j,k=0 j,k=0 Z=Max(j+l,fc-|-l) 


e-1 


Min(e, j+k+1) 

'Y. o,j+k-i+i h . (A9) 

j,k=0 /=Max(j-l-l,/i;+l) 


Y 


For polynomials of the same degree {d = e) the result simplifies to directly yield 

Min(d, _ 7 +fc+l) 

Bjkif,g)= Y {(^ibj+k-i+i - aj+k-i+ik). (AlO) 

/=M ax (_ 7 '+1 ^/c+1) 

Assuming the more general d > e case, the resultant is computed as the determinant of the 
Bezout matrix multiplied by a non-zero factor: 


Res,(/,j) = (-l)‘*<‘‘-‘>/2a5-idetB(/,9). (All) 


Hence, the condition for a common root from equation (A5) is equivalent to det!?(/, g) = 0. 


The Bezout matrix has the advantage of being smaller, while the Sylvester matrix is 
usually sparser. If any particular symbolic-manipulation software package fails to compute 
the determinant of one of them, one can try using the other, before resorting to numerical 
computations. 
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B. Spurious results of the resultant method 


The polynomial pres{w) defined by eqnation ( pT| has a set of roots wi, ... Wm- The 
polynomial Pconj{w) defined by eqnation (12) then has roots The search for a 


common nnit root nsing eqnation (13) assnmes that for at least one of the roots Wj = w 


-1 


and hence it lies on the nnit circle. However, eqnation (13) is satisfied even by a non-nnit 
root Wj if it is eqnal to the inverse conjngate of a different root, i.e., Wj = for some other 


Wk 7 ^ Wj. This implies that Wk is another non-nnit root satisfying eqnation (13). In case p. 


has only real coefficients (e.g., in the case of the two-point-mass lens or the triple-lens models 
from 


5.2 


5.5) and Wj lies off the real axis, Wj and Wk are two more non-nnit roots of 


Pres- The total nnmber of such spurious solutions depends on the degree of the polynomial 
Pres, which is < 3n — 3, as discussed in 


3.2 


In the case of the two-point-mass lens, the polynomial given by equation ( |l7| is of third 
degree and has real coefficients. Spurious solutions off the real axis would require Pres to 
have at least four different roots {tci, W 2 = Wi^, wi, W 2 = hence they do not occur. 

A spurious solution wi along the real axis would require pres to have a set of real roots 
{tci, W 2 = Wi^, Wa} and Pconj a corresponding set of roots could then 

write 

Presiw) = {w- Wi){w - W^^){w - W3) . (Bl) 

Expanding this expression and comparing it term-by-term with equation dl?] ), we find that 
the only possible real solution is wi = W 2 = 1 , which lies on the unit circle and thus is 
not spurious. Moreover, this double-root solution requires either /r = 0 or /r = 1, and 
thus it corresponds to the single-lens limit rather than a genuine binary. The nonexistence 
of real spurious solutions can also be seen from the discriminant of Pres computed from 
equation (ED' The result is always negative for /i 7 ^ 0, 0.5, 1. Therefore, except for these 
values Pres has one root on and two roots off the real axis. For the only other non-degenerate 
case, /i = 0.5, the roots are {—2s^, —2s^, s^/4}. We see that even in this case there is 
no reciprocal pair of roots and hence no spurious solution. For the two-point-mass lens all 


solutions of equation (13) thus correspond to topology changes of the critical curve. 


In all the triple-lens cases analyzed in § 5^ - § |5.5, spurious solutions do occur. In the 


LS model the 6 th degree saddle-point and critical-curve polynomials in fact can be studied 
as 3rd degree polynomials in z^. The first resultant then leads to a 3rd degree Pres{w), which 
can be shown to have a pair of real spurious roots wi = w^^ when the central mass p < 0 . 2 . 
In the three other models pres{w) is of 6 th degree. Regions in parameter space with spurious 
solutions can be found for each of them. 
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C. Cusp and morph curve polynomials of the LS model 


We present here the explicit form of the cusp and morph curve polynomials for the 
simplest triple-lens model studied in this work, the linear symmetric model from 


5.2 


The cusp curve is generally dehned by equation (27), with n = 3 for the triple lens. 
We hrst multiply the equation by {z — zi)^(z — Z 2 )^(z — ^ 3 )®. Next we substitute the posi¬ 
tions and masses of the LS-model components, {zi,Z 2 ,Z 3 } = {—s,0,s} and {jui, jU 2 , jus} = 
{(1 — /i)/2,/i, (1 — /i)/2}. We obtain the parametric polynomial form of the cusp curve. 


Pcusp(^^') 

(1 - A) + 3 [2 - 4/i + A(3/r - 1) ] s^z^° + 3 [3 - 10/i + 12 / 1 ^ - A(1 - 5/r + V) ] sS® 
+ [4/i(4 - 9/r) - A(1 - 3/i + V - 27/^3) ] - 3 [2 - 7/r + A(1 - 5/r + V) ] s^z^ 

—3 [2 -|- A(1 — 3/i) ] + (1 ~ = 0 , (Cl) 

where A > 0 is a real non-negative parameter. 


In a similar manner we convert the morph-curve equation (28) for n = 3, using exactly 
the same steps. After multiplication by (z — zi)^(z — Z 2 )^(z — zs)^ and substitution of the LS 
model positions and masses, we obtain the parametric polynomial form of the morph curve. 


Pmorphi^^^ 


-i Tz^^ + [ 1 - /i - 6 i r(l - 2/i) ] 5^^^° - [2 - 8/i + b/i^ + 3 i r(3 - 10/i + ] sV 

-|-[ 1 — 12/i -|- 11/i^ — 4 ir/i(4 — 9/i) ] s^z^ + [4 — 4/i-|-3ir(2 — 7/i)]/i s^z'^ 

+ [1 -/i + 6ir/i]/is^V -iT/i^s^^ = 0, (C2) 


where T is a real parameter. 

We see that in the LS model both polynomials are of sixth degree in z'^. Intersections of 
the cusp curve given by equation (Cl) and the critical curve given by equation (33) identify 
cusp images along the critical curve. Intersections of the morph curve given by equation (C2) 
and the cusp curve identify the images of caustic-metamorphosis points in the image plane, 
as discussed in § |^ 
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Fig. 1.— Configurations of the analyzed triple-lens models. Lens components are marked by 
numbered circles and labelled by their masses; black circles indicate equal mass. Model pa¬ 
rameters: Linear Symmetric - central mass fi, separation s of neighbors; Linear Asymmetric 
- relative position p of central component, half-separation s of outer components; Triangular 
Equilateral - mass p of component 1, side length s; Triangular Isosceles - vertex angle 6, 
leg length s. 
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Fig. 2.— LS model parameter-space division. Left panel: by critical-curve topology (marked 
by letters), with black boundary given by equation (36), orange by equation (37), cyan by 
Right panel: by total number of cusps on caustic [8 cusps - gray; 12 - red; 
brown]. Region near lower left corner of right panel is blown up in Figure]^ 


equation (38). 
16 - cyan; 20 
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Fig. 3.— LS model: detail of parameter-space division by total cusp number from Figure 
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Fig. 4.— LS model: gallery of topologies of critical curves (blue) and corresponding caustics 
(red) with lens positions marked by black crosses. Letters correspond to regions marked 
in Figure Caustic subregions and lens parameters [p, s] of examples: Ai [0.1,1-5], B 2 
[0.2, 0.5], Cl [0.04,0.5], D 2 [0.2, 0.4], E 3 [1/9,0.35], Fi [0.07,0.35]. 
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Fig. 5.— LS model: gallery of critical curves (blue) and corresponding caustics (red) in 
transitions between the regions marked in FigureLetters denote all adjacent regions. Lens 
parameters [/i, s] of examples: AB [1/3,1-6777], BC [0.05,0.53075], BD [1/3,0.54593], BE 
[0.1, 0.38852], CF [0.05, 0.36951], DE [0.15, 0.33], EF [0.08, 0.34866], BCEF [0.07530, 0.37866], 
BDE [1/5, 5-1/2], 
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Fig. 6 .— LS model: gallery of caustics with blown-up details from 13 subregions of right 
panel of Figure plus two transition caustics. Labels mark cusp-number subregions, with 
numbers assigned within given region from top to bottom and left to right in Figure 
Lens parameters [/i, s] of examples: Ai [0.1,1.5], A 2 [0.95,1.45], Bi [0.2, 0.8], B 2 [0.2, 0.5], 
B 3 [0.04,0.6], Cl [0.04,0.47], Di [0.7,0.65], D 2 [0.2, 0.4], Ei [0.15,0.36], E 2 [0.155,0.35], E 3 
[1/9,0.35], E 4 [0.08,0.35], Fi [0.07,0.36]; BCEF [0.07530,0.37866], BDE [l/ 5 , 5 -i/ 2 ]_ 
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Fig. 7.— LS model: Jacobian contour plots for BCEF boundary intersection at /i ~ 0.0753 
(left column), for double maxima at p = 1/9 (central column), and for BDE monkey saddles 
at /i = 1/5 (right column). Lower row includes cusp curves (orange) and morph curves 
(green). Contour values differ from column to column. Bold black contour: critical curve 
for s = 0.5; crosses: lens components; pluses: Jacobian saddle points. 
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Fig. 8.— LA model parameter-space division. Left panel: by critical-curve topology (marked 
by letters), with black boundary given by PresiX) = 0, cyan by 15th degree polynomial in 
(see § 5.3). Right panel: by total number of cusps on caustic [8 cusps - gray; 10 - blue; 12 
- red; 14 - green; 16 - cyan; 20 - brown]. 
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Fig. 9.— LA model: gallery of topologies of critical curves (blue) and corresponding caustics 
(red) with lens positions marked by black crosses. Letters correspond to regions marked in 
Figurej^ Caustic subregions and lens parameters [p, s] of examples: Ai [0.2,1.2], Bi [0.4,1.5], 
Cl [0.5,1.725], Ds [0.5,0.75], Ei [0.7,0.75], Fa [0.5,0.45]. 











Im[z] Im[z] 


57 



Re[z] 


Re[z] 




-1 0 1 

Re[z] 


Fig. 10.— LA model: Jacobian contour plots in regime of binary lens with distant companion 
at p = 0.18 (left column), for the nearly symmetric p = 0.45 (central column), and the 
symmetric case p = 0.5 (right column). Lower row includes cusp curves (orange) and morph 
curves (green). Contour values differ from column to column. Notation as in Figure]^ 
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Fig. 11.— TE model parameter-space division. Left panel; by critical-cnrve topology 
(marked by letters), with black boundary given by PresiX) = 0, orange by Pres{—X = 0^ 
cyan by 3rd degree and green by 12th degree polynomials in (see § 5. 
total number of cusps on caustic [ 10 cusps - blue; 12 - red; 14 - green; 16 
yellow; 20 - brown]. Region near lower right corner of right panel is blown up in Figure 12 


Right panel: by 
cyan; 18 - 
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Fig. 12.— TE model; parameter-space division by total cusp number. Left panel: detail 
from Figure Right panel: detail from left panel near close limit of double-maximum 


configuration. 
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Fig. 13.— TE model: gallery of topologies of critical curves (blue) and corresponding caustics 
(red) with lens positions marked by black crosses. Letters correspond to regions marked in 
Figure 11, Caustic subregions and lens parameters [/i, s] of examples: A 4 [1/3,1-55], B 3 


[0.15,1.3], Cl [0.7685,1.0], Di [1/3,1.2], E 2 [0.7685,0.68], F 2 [0.9,0.68], Gs [1/3,0.65], Hi 
[0.15,0.68], Ii [0.6,0.68], J 3 [0.7685,0.64]. 
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Fig. 14.— TE model: Jacobian contour plots for equal masses /i = 1/3 (left column), for 
double maximum at /r ~ 0.768 (central column), and for monkey saddle at fi = 8/9 (right 
column). Lower row includes cusp curves (orange) and morph curves (green). Contour values 
differ from column to column. Insets in lower central and right panels: zoomed-out plots of 
critical, cusp, and morph curves. Notation as in Figure 
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Fig. 15.— TI model parameter-space division. Left panel: by critical-curve topology (marked 
by letters), with black boundary given by Pres(l) = 0, orange by presi,—^) = 0, cyan by 15th 
degree polynomial in (see § 5.5). Right panel: by total number of cusps on caustic [8 
cusps - gray; 10 - blue; 12 - red; 14 - green; 16 - cyan; 18 - yellow; 20 - brown]. Regions 
near upper left and lower right corners of right panel are blown up in Figure 
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Fig. 16.— TI model; parameter-space division by total cusp number. Left panel: detail 


from upper left part of Figure 15 Right panel: detail of lower right corner of Figure 15 
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Fig. 17.— TI model; gallery of topologies of critical curves (blue) and corresponding caustics 
(red) with lens positions marked by black crosses. Letters correspond to regions marked 
in Figure 15 Caustic subregions and lens parameters [0,5] of examples: A 2 [0.05 tt, 2 ], 


B 2 [0.l7r,2], C 2 [0.27r,2], Dg [7r/3,1.55], E 2 [0.l7r,1.2], Fi [0.27r,1.2], G 2 [0.77r,1.2], H 2 
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Fig. 18.— TI model: Jacobian contour plots for off-axis butterflies at 6 ^ 0.220 tt (left 
column), for Chang-Refsdal violating 6 = 7r/2 (central column), and for double maximum 
at 6^ ~ 0.708TT (right column). Lower row includes cusp curves (orange) and morph curves 
(green). Contour values differ from column to column. Insets in lower panels: zoomed-out 
plots of critical, cusp, and morph curves. Notation as in Figure 


































































Table 1. 


Critical-curve Topology Occurrence in Triple-lens Models 


Model 

o 

oo 

o 

ooo 

OO 

© 

©o 

© 

© 

Linear Symmetric (LS) 

B 

- 

- 

A 

- 

C,E 

- 

- 

D,F 

Linear Asymmetric (LA) 

D 

B 

- 

c 

- 

E 

A 

- 

F 

Triangular Equilateral (TE) 

C 

B 

D 

A 

- 

E,H 

- 

I,J 

F,G 

Triangular Isosceles (TI) 

G 

C 

F 

D 

B 

E,J 

A 

I,L 

H,K 


Note. — Topology regions of the model parameter spaces are marked by letters, as defined in 
Figure for the LS model, Figure for the LA model. Figure [TT] for the TE model, and Figure [T^ 
for the TI model. 
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